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Chapitre 1 

Introduction aux équations de 
Lagrange 


1.1 Equations de Lagrange pour une particule 

1.1.1 Equations de Lagrange 

Considérons le cas particulier d’une particule astreinte à se déplacer, sans frottement, sur 
une courbe plane contenue dans le plan xOy. La courbe sur laquelle est astreinte se déplacer la 
particule de masse m, est le lieu des points dont les coordonnées vérifient les relation : 


z = 0 

f{x,y) = 0 


La première relation correspond au plan xOy . La seconde relation représente l’équation de 
la trajectoire dans ce plan. Ces deux relations définissent les équations des liaisons appelées 
souvent liaisons. Le nombre de degrés de liberté est égal au nombre de coordonnées qui repré- 
sentent la position de m ( trois dans le cas général) moins le nombre de liaisons ( deux dans le 
cas présent ). La particule possède donc un degré de liberté. Il faut choisir une variable q pour 
repérer sa position. Cette variable est appelée coordonnée généralisée. Il est possible d’exprimer 
le vecteur position f de la particule en fonction de la coordonnée généralisée q pr la relation : 
r = f(q). 

Soit F la résultante de toutes les forces agissant sur la particule. La relation fondamentale 
de la dynamique s’écrit : 


cPr dv 

r = m—rr = m— 

dt 2 dt 


où v = — est la vitesse de la particule. 

Soit 8W le travail fourni par la force F lors d’un déplacement infinitésimal 5r : 


ÔW = F -ôf 
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Le déplacement Sr peut s’écrire en fonction de la variation Sq de la coordonnée généralisée 


Q ■ 

dr 

Sf= —Sq 
dq 

Dans ce cas le travail SW peut se mettre la forme : 

-> dr 

SW = F ■ —Sq 
oq 

On appelle force généralisée conjuguée de q, ou ç-composante de la force, la quantité F q 
définie par : 

SW p> dr 
q = li = "dq 

Par conséquent SW s’écrit : 

SW = F q Sq 

En tenant compte de la relation fondamentale de la dynamique, cette expression peut éga- 
lement s’écrire : 


dv dr 

SW = m— ■ - —Sq 
dt dq 


D’autre part 


Sachant que 


d 

^ dr 

dv dr ^ d 

dr 

dt 

V ' dq. 

dt dq dt 

dq 


d 

dr 

d 

dr 

dt 

dq 

dq 

dt 


dv 

dq 


on obtient 

dv dr d 

dt dq dt 

Le vecteur vitesse v. peut aussi s’écrire : 


v ■ 


dr 

dq 


— v ■ 


dv 

dq 


D’où la relation : 


et 


^ dr drdq dr . 

dt dq dt dq ^ 


dr dv 
dq dq 


dv dr d 

dt dq dt 


v ■ 


dv 

dq 


— v ■ 


dv 

dq 
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1.1 Equations de Lagrange pour une particule 


Sachant que 


et que 


on obtient 


d \ 1 2 ~| <9 H 1 _ dv 

d~q 2 V = OÏi 2 V ' V =V 'd7i 


d f 1 2 1 d n 1 _ dv 

dï 2 V = ~dq 2 V ' V =V '~dq 


clv dr d \ d \1 


dt dq dt | dq |_2 ? 
L’expression du travail SW peut alors s’écrire : 


d fl 


dq 2° 


f d \ d \ 1 2 1 1 d \ 1 2 1 1 

Sw = m U 94 2” \\-¥ q [2 V I Sq 


Si on note T = ^mv 2 l’énergie cinétique de la masse m , on obtient finalement : 

ÂU .jd\dT] dT 1 
' { dt [ dq \ dq J q 

On obtient finalement les deux expressions équivalentes du travail SW 

( d \dT] dT\ 

\dt[aï\--a q r 9 = FA 

On en déduit l’équat-ion de Lagrange pour un système à un degré de liberté : 

d \dT]dT_ 
dt dq dq q 


1.1.2 Cas des systèmes conservatifs 

Dans les systèmes conservatifs, la force appliquée au système dérive d’un potentiel U et elle 
s’écrit : 

F = J_E_ 

q dq 

L’équation de Lagrange devient alors : 

d fdTl dT dU 


dt [dq \ dq dq 


Généralement l’énergie potentielle U ne dépend pas de la vitesse, c’est-dire que — — = 0. 

dq 

L’équation de Lagrange peut alors s’écrire : 

d_ \ d(T-U) 1 _ d (T - U) _ 
dt dq dq 
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On introduit la fonction de Lagrange ( ou lagrangien du système ) qui est la différence de 
l’énergie cinétique et de l’énergie potentielle : 

C = T — U 

D’où la forme de l’équation de Lagrange dans le cas d’un système conservatif : 

d_ raci _ oc _ 

dt dq dq 


1.1.3 Cas des forces de frottement dépendant de la vitesse 

Equation de Lagrange 

Considérons une situation physique dans laquelle la particule est soumise à des forces de 
frottement de viscosité dont la résultante / est de la forme : 

f = -av 


Pour calculer la force généralisée f q correspondante, nous utilisons la définition du para- 
graphe précédent : 


U 



= —a 


dr 2 dq 
dq dt 


Cette dernière expression peut se mettre sous la forme : 



cl \dc 1 _ de 

dt dq dq 
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Fonction dissipation 

Calculons le travail ÔWf fourni par la force de frottement pendant un intervalle de temps 
St pour un déplacement Sf : 

SWf = f ■ Sf = — a v 2 St 

La quantité de chaleur SC) gagnée par le système en interaction avec la particule, est telle 
que : 


SQ = a v 2 St 


SQ 


Soit P d = la puissance dissipée par les forces de frottement sous forme de chaleur : 
St 

P d = a v 2 

Cette puissance dissipée peut être exprimée en fonction de q, par : 

2 


Pd = Ot 


dr 

2 

drdq 

dt 

= Oi 

dq dt 


= (3 q 


•2 


Par définition, la fonction dissipation est égale à la demi-puissance dissipée : 
La ^/-composante f q de la force de frottement peut alors s’écrire : 


fç = - 


do 

dq 


L’équation de Lagrange s’écrit alors : 

d \dC 
dt dq 


dC dv_ 
dq ^ dq 


1.1.4 Cas d’une force extérieure dépendant du temps 

Considérons le cas plus général d’un force extérieure dépendant du temps agissant sur un 
système qui est le siège de forces de frottement qui dérivent d’une fonction dissipation V. Soit 
F eq la (/-composante de la force extérieure. Dans ce cas l’équation de Lagrange peut s’écrire 
sous l’une des deux formes équivalentes suivantes : 


d 

~dc 

de 

dt 

dq 

dq 

d 

' dC 

dC 

dt 

dq 

dq 


F eq ~ P q 


dv 

dq 


P 


e,q 
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1.2 Système à plusieurs degrés de liberté 


Dans le cas général d’un système à plusieurs degrés de liberté, il y a autant d’équations 
de Lagrange que de degrés de liberté. Ainsi, si le système possède N degrés de liberté, il est 
nécessaire d’avoir N coordonnées généralisées Qi ( i = 1, 2, N) ; nous aurons ainsi N équations 
de Lagrange : 

d 
dt 


dC dC dV ^ .. 1 n 

Mi (î “ 1 ’ 2 ’-"^ ) 


La q— composante de la force généralisée extérieure est définie par : 


F = 

1 e,q i 


SW 


Sqi 


Sqi^O 


Dans cette expression SW représente le travail des forces extérieures résultant d’une variation 
Sqi de la coordonnée q t telle que les coordonnées soient constantes {Sq -^.-, = 0). 
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Chapitre 2 

Oscillations libres des systèmes à un 
degré de liberté 


2.1 Oscillations non amorties 

2.1.1 Oscillateur linéaire 

Un système oscillant à un degré de liberté est habituellement repéré à l’aide d’une coor- 
donnée généralisée g qui est l’écart par rapport à la position d’équilibre stable. Le mouvement 
vibratoire est dit linéaire s’il est régi par une équation différentielle harmonique de la forme : 

q + UJ 2 0 q = 0 

Cette équation est appelée équation différentielle de l’oscillateur harmonique simple. 


2.1.2 Energie Cinétique 

Dans le cas d’un système à un degré de liberté, constitué d’une masse m dont la position 
est répérée par la coordonnée généralisée g, l’énergie cinétique s’écrit : 


1 2 1 

T = -mv = -m 
2 2 


df 

dt 


1 

- -m 


dfdq 
dq dt 


1 

- -m 


df 

dq 


L’énergie cinétique d’un système à un degré de liberté est fonction de q et q . Elle peut 
s’écrire sous la forme : 


T = \a(q) q 2 

où a(q) est une fonction de la coordonnée généralisée g, définie dans le cas étudié par : 


a (g) = m 
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En faisant un développement limité de a(q) au second ordre en g, au voisinage de q = 0 , on 
obtient : 

.. 1 ... da 1 d 2 a 2 , 2 

T(q ' q)= 2[ a{0)+ T iq J + 2Së ,J + -\ q 

En limitant l’approximation au second ordre, on obtient : 


rri • 2 

T = -a 0 q 

où «o est une constante égale à a (0) . 

2.1.3 Energie potentielle 

Les oscillations se font autour de la position d’équilibre stable q = 0 caractérisée par 
dU 

——(q = 0) = 0. Il est toujours possible , lorsque les écarts par rapport à la position d’équilibre 
dq 

sont faibles, de faire un développement en série de Taylor de U (q) au voisinage de la position 
d’équilibre q = 0. En négligeant les puissances de q d’ordre supérieur à deux, on obtient : 


U(q) = U( 0) 


dU 1 d 2 U 2 


q = 0 correspond à un minimum de U ( q ) pour lequel 

dU d 2 U 

— = 0 et —~y > 0 

d( l q =o d( l q =o 

Si on choisit l’origine de l’énergie potentielle à cette position d’équilibre (1/(0) = 0) , l’énergie 
potentielle U ( q ) peut s’écrire sous une forme quadratique : 

U(q) ^ ^ 2 


avec : b n = 


dq 2 „=n 


2.1.4 Equation différentielle 


L’équation de Lagrange s’écrit : 


d p£l _ dC _ 
clt [dq \ dq 


Ce qui permet dobtenir l’équat-ion différentielle de l’oscillateur harmonique simple où : 


d 2 U 

dQ 2 n —r 


UJ 0 = — = 
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Les oscillations d’un système vibratoire s’effectuent autour d’une position d’équilibre stable. 
Pour des oscillations de faible amplitude autour de la position d’équilibre, tous les mouvements 
vibratoires peuvent être assimilés à des vibrations linéaires et l’énergie potentielle peut alors 
être approximée par une forme quadratique de la coordonnée q, tandis que l’énergie cinétique 
peut être approximée par une forme quadratique en q. 

2.1.5 Résolution de l’équation différentielle de l’oscillateur harmo- 
nique simple 

L’équation différentielle de l’oscillateur harmonique simple s’écrit : 

q + u 2 0 q = 0 

La solution d’une telle équation est une fonction sinusoïdale du temps 

q(t) = A cos (cu 0 1 + <p) 

où A représente l’amplitude des oscillations, (p est la phase initiale. 

Il est important de remarquer que la pulsation propre oüq ne dépend que des éléments qui 
constituent le système physique étudié (masse, ressort, etc...) tandis que l’amplitude A et la 
phase initiale ip sont calculées à partir des conditions initiales : 

f q(t = 0 ) = q 0 
1 q(t = 0 ) = q 0 

Enfin l’amplitude des oscillations d’un oscillateur harmonique libre ne dépend pas du temps. 
De telles oscillations sont dites non amorties. 

Il faut néanmoins remarquer qu’au delà d’une certaine amplitude la vibration devient non 
linéaire. Il s’ensuit d’abord une modification de la période des oscillations et ensuite un chan- 
gement de la nature du mouvement. 

2.2 Oscillations libres des systèmes amortis à un degré 
de liberté 

Dans le paragraphe précédent, nous n’avons pas tenu compte de certaines réalités physiques. 
En effet, nous n’avons pas pris en compte les forces de frottement qui sont à l’origine de la perte 
d’énergie mécanique du système sous forme de chaleur. Dans ce paragraphe, nous allons tenir 
compte de ces réalités en nous limitant- toutefois au cas simple où les pertes sont dues à des 
frottements visqueux pour lesquels les forces de frottement, qui s’opposent au mouvement, sont 
proportionnelles à la vitesse . 
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2.2.1 Equation de Lagrange pour les systèmes dissipât ifs 

Rappelons l’équation de Lagrange associée à un système à un degré de liberté dont l’évo- 
lution au cours du temps se ramène à l’étude de la coordonnée généralisée q 

d_ 

clt 

Fq représente la composante suivant q de la résultante des forces généralisées qui ne dérivent 
pas d’un potentiel. 

Nous nous intéressons au cas particulier des forces de frottement définies par la force géné- 
ralisée 


dC _ dC _ 
dq dq q 


F q = fg = -P4 

où 3 est une constante réelle positive. 

L’équation de Lagrange s’écrit alors dans ce cas : 

d dC 
dt dq 



2.2.2 Cas particulier des oscillations de faible amplitude 

Nous avons montré dans le chapitre précédent que dans ce cas, la fonction de Lagrange 
s’écrit sous la forme : 

r 1 -2 2 

£ = 2 aq ~ 2 bq 

L’équation différentielle du mouvement s’écrit alors : 

a q + bq = ~3q 


C’est une équation différentielle du second ordre à coefficients constants qui peut se mettre 
sous la forme : 

q + 28q + ulq = 0 

où S est un coefficient positif, appelé facteur (ou coefficient) d’amortissement et défini par : 


lü(j est la pulsation propre définie par 


S =7T- 
2 oq 
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2.2.3 Résolution de l’équation différentielle 

La solution de l’équation différentielle dépend de la valeur de S par rapport à cu 0 : 

- Si S > uj() . on dit que le système est suramorti ou apériodique. 

- Si 5 = ujq , on dit que l’on a un amortissement critique. 

- Si 5 < ujo , on dit que le système est sous-amorti ou pseudopériodique. 

Cas où le système est suramorti (S > ujo) 

La solution de l’équation différentielle s’écrit dans ce cas : 

g(t) = A 1 e[^^ f + d 2e H + ^‘ 

Ai et A 2 sont des constantes d’intégration définies par les conditions initiales. La figure 
ci-dessous représente g en fonction du temps dans le cas particulier où q( 0) = go et g(0) = 0. 
q(t) est une fonction qui tend exponentiellement (sans oscillation) vers zéro. 



Régime fortement amorti : variation de q en fonction du temps 


Cas de l’amortissement critique (ô = cu 0 ) 

La solution générale de l’équation différentielle est de la forme : 

q(t) = ( Ai + A 2 1) e 5t 

Dans le cas particulier où ç(0) = go et g(0) = 0, 

q(t) = g 0 (1 + St) e~ 5t 

q(t ) est encore une fonction qui tend vers zéro sans oscillation lorsque le temps augmente. 
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Amortissement critique : variation de q en fonction du temps 


Cas où le système est sous-amorti (5 < o; 0 ) 

La solution générale de l’équation différentielle est de la forme : 

q(t) = A e~ St cos ((jüa t + 4>) 


avec üj,\ = a/cUq — ô 2 ; A et 0 sont deux constantes d’intégration déterminées à partir des 
conditions initiales. Dans le cas particulier où q( 0) = qo et q(0) = 0, on obtient : 


A 


UJ 0 

1 

UJ A 


-qo 


arctan — 



Système faiblement amorti : variation de q en fonction du temps 
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La figure représente les variations de q(t) au cours du temps. On remarque que q(t) est 
enveloppée par les deux fonctions exponentielles : q 0 e~ St . 

Le lieu des maxima est obtenu en résolvant q(t) = 0. Ce qui donne : 

tan(cüAt + 0) = = tan (0) 

UJa 

d’où l’on tire l’instant- t n correspondant au n-ième maximum : 

2vr 

t n — Tl 

UJA 


Les maxima de q(t) sont séparés par des intervalles réguliers égaux à 



UJ A 


Ta est appelée la pseudo-période. On remarque que , en plus de la diminution de l’amplitude 
des oscillations au cours du temps, l’un des effets de l’amortissement est l’augment-ation de la 
période des oscillations. Pour des systèmes faiblement amortis (5 « uj 0 ), on peut remarquer 

271 

que uja — <^o et- que la pseudo période est peu différente de la période propre : Ta — T 0 = — . 

Cd 0 


2.2.4 Exemples 

Système mécanique en translation 

Amortisseur mécanique Un amortisseur mécanique est constitué d’un élément mobile à 
l’intérieur d’un récipient contenant- un fluide visqueux. 


a 






l> 





Amortisseur 


La force de frottement- / agissant sur la partie mobile repérée par aq, est donnée par 

fx = -a (ii - ;r 2 ) 

où (ii — i 2 ) représente la vitesse relative des deux éléments qui constituent l’amortisseur. 
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Equation différentielle du mouvement Considérons le cas d’une masse m oscillant ver- 
ticalement et reliée à un bâti fixe par un ressort de raideur k et un amortisseur de coefficient 
de frottement visqueux a. Repérons par x l’écart de la masse m par rapport à la position 
d’équilibre. L’équation de Lagrange s’écrit : 

d_ ïdc 

clt dq 

Sachant que T = \mx 2 et que U = \kx 2 , la fonction de Lagrange se met sous la forme 

r 1 • 2 1 7 2 

L = -mx — -kx 
2 2 

On obtient l’équation différentielle du mouvement 

mx + ax + kx = 0 

Cette équation peut s’écrire sous la forme d’une équation différentielle du second ordre à 
coefficients constants 


dC 

dq 


= —ax 


x + 25 x + lü^x = 0 

où cüq et 5 sont des constantes positives appelées respectivement la pulsation propre et le 
facteur d’amortissement ; elles sont données par les relations suivantes 


5 


a 

2 m 


et cuq 



Système mécanique en rotation 

Considérons un pendule simple constitué par une masse ponctuelle m reliée à un axe de 
rotation fixe par une tige de masse négligeable, qui effectue des oscillations de faible amplitude 
dans un plan vertical. Ce dispositif est un système à un degré de liberté dont on repère la 
position par l’angle 0 par rapport à la verticale. Le système est immergé dans un fluide dont les 
forces de viscosité se ramènent à une force s’exerçant sur la masse m et donnée par la relation 

/ = -otv 


L’équation de Lagrange qui régit le mouvement d’un tel dispositif s’écrit : 

d 
dt 

L’équation différentielle du mouvement s’écrit alors : 

ml 2 9 + al 2 9 + m g 1 6 = 0 


OC 

dq 


dC 

dq 


= —al 2 9 
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9 + 269 + uj 2 0 9 = 0 

où la pulsation propre co 0 et le facteur d’amortissement S sont respectivement donnés par : 


UJ 0 



et 


5 


a 

2 m 
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Chapitre 3 

Oscillations forcées des systèmes à un 
degré de liberté 


3.1 Equation différentielle 

Rappelons la forme générale de l’équation de Lagrange pour les systèmes à un degré de 
liberté : 


d \ÔC 
dt dq 


dC dV _ p 

, , h r-, . -C qext 

oq oq 


où F qext est la force généralisée associée à F ext et où la fonction dissipation est V = 

2^ 2 ' 

Pour les oscillations de faible amplitude, la fonction de Lagrange pouvait se mettre sous une 
forme quadratique de q et q 

r 1 - 2 2 

£ = 2 °° d _ 2° q 

D’où l’équation différentielle du mouvement 


°o q + (3 q + b 0 q — F qext 

Cette équation peut se mettre sous la forme d’une équation différentielle du second ordre à 
coefficients constants, avec second membre 

q + 2Sq + ulq = A(t ) 


avec 


5 


A 

2oq ’ 


W 0 



et A(t) 


Fqext 

«O 
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3.2 Système masse-ressort-amortisseur 



Système masse-ressort-amortisseur 


Considérons l’exemple mécanique de la figure ci-dessus soumis à une force extérieure F(t) 
appliquée à la masse m. Calculons la force généralisée F x conjuguée de la coordonnée x. Pour 
cela nous pouvons utiliser l’une des deux méthodes suivantes : 

- Soit calculer le travail dW de la force F(t) pour une variation dr de son point d’application 

dW = F ■ dr = F dx 


On en déduit la ^-composante de la force extérieure 

f = FW 

- Soit utiliser la définition de la force généralisée 

— * f)r 

F . = F -is = m 

L’équation différentielle du mouvement s’écrit alors 

x + 2ôx + uj^x = A(t ) 


avec : 


S = a/2m , ujo = \j - — et A(t) = F(t)/m 
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3.3 Solution de l’équation différentielle 

La solution de cette équation différentielle du second ordre est égale à la somme de la solution 
de l’équation sans second membre (ou solution homogène) x h (t) et d’une solution particulière 
de l’équat-ion avec second membre Xp(t) : 

x(t) = x H (t ) + x P (t) 

Nous avons déjà étudié l’équat-ion sans second membre xp(t) et nous savons que cette 
solution contient dans tous les cas le terme exponentiel e 5t . Après un temps t supérieur à 3 /ci ou 
4/à, le terme e 5t devient très petit et la solution homogène est alors pratiquement nulle. Il ne 
subsistera que la solution particulière de l’équation avec second membre. L’intervalle de temps 
pendant lequel la solution homogène est non négligeable est appelé le régime transitoire. A la 
fin de ce régime transitoire commence l’intervalle de temps pour lequel la solution homogène 
est quasi- nulle et pour lequel la solution x(t) ~ x p (t) ; ce régime est appelé régime permanent 
ou stationnaire. 


3.3.1 Cas particulier où A(t) = Aq cos(12t) 

a) Calcul de la solution permanente l’aide de la méthode des nombres 
complexes 

Pour t suffisamment grand, nous pouvons considérer que la solution transitoire s’est- annulée 
et que la solution x(t) s’identifie alors avec la solution particulière : x(t ) ~ Xp(t). Par commodité 
de notation l’indice p est sous-entendu dans ce qui suit. La méthode des nombres complexes 
permet de calculer aisément la solution stationnaire. 

Soit le déplacement complexe représenté par le nombre complexe X = avec X = 

X 0 e J:p . Nous pouvons considérer, en outre, que A(t) = A 0 cos (Qt) constitue la partie réelle du 
nombre complexe A = Aq e j ' H . L’équation différentielle se transforme en une simple équation 
algébrique en fonction de l’amplit-ude complexe N : 

[(wg -fi 2 ) + j2SQ]X = Ao 

dont la solution est : 


— (u$-n 2 )+j2sn 

D’où l’on tire l’amplit-ude X 0 et la phase <p : 


\/(v$-n 2 f + 4ô 2 n 2 

2<5D 

<p = — arctan — » — 

CUq — il 2 
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b) Etude des variations de l’amplitude et de la phase en fonction de 
la pulsation de l’excitation 

dX 

Le maximum de l’amplitude est obtenu pour la valeur de fi qui annule — — . 

ail 

Il existe un maximum à la pulsation fi R = a/cUq — 2 8 2 seulement si l’amortissement est 
suffisamment faible pour que S < ujq/\/2. A cette pulsation appelée pulsation de résonance, on 
dit que le système entre en résonance et l’amplitude X 0 est maximale ; elle vaut : 


X, 


Ar 


0 max 


2 ôy/ul ~ Ô 2 

La figure représentant les variations de X 0 en fonction de la pulsation d’excitation fi est 
appelée courbe de résonance en amplitude. On remarque qu’à la pulsation uj 0 , le déphasage (p 


7T 


est égal à — — , et qu’à la résonance p = — arctan 

A 


\Ju%- 2 ô 2 
S 



c) Etude de la résonance pour les faibles amortissements 

Dans le cas des faibles amortissements ( 8 « oj 0 ), la fréquence de résonance est très peu 
différente de la pulsation propre, fi/,> ~ co 0 . Dans ce cas, l’amplitude de vibration à la résonance 
Ao max est égale à : 

y - 

Av. 0 max n e 

ZôüJ o 

Pour les faibles amortissements, X (} rnax est donc inversement proportionnel à 8. 


d) Etude de la vitesse 

En notation complexe, la vitesse s’écrit : 

V{t) = ^= ;jÜX = Xe jQt 
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où l’amplitude complexe de la vitesse est définie par 


X = jVtX 


jttA o 

{ujl-Q?)+j25Ü 


L’étude des variations de l’amplitude de la vitesse en fonction de la pulsation d’excitation 
montre que, quelle que soit la valeur de <5, la résonance en vitesse est obtenue pour Q = u (J 
(voir figure ci-dessous). La valeur maximale de l’amplitude de la vitesse vaut dans ce cas : 


di-max — X(CJ 0 ) 


Ao 
2 ~ô 



n/2 


0 


-n/2 



Courbe de résonance de la vitesse Déphasage ip de la vitesse en fonction de Q 


e) Bilan énergétique 

Soit Pp(t) la puissance instantanée fournie par la force extérieure F(t) au système. En 
régime permanent, on obtient : 

Pp(t) = F(t) x(t ) = F 0 X 0 cos (Qt) cos (Qt + ip) 

Soit < Pp > la valeur moyenne sur une période de Pi-(t) : 

< P f >= -/ 0 X 0 cos (VO 

En tenant compte de l’expression de X 0 en fonction de F 0 , on obtient : 

1 • 9 

< Pf >= 2 aX o 

Comparons cette valeur à la valeur moyenne < Pp > de la puissance dissipée par les forces 
de frottement de viscosité. La valeur instantanée de cette puissance dissipée s’écrit : 

Po(t) = ax 2 = aXl cos 2 (Vit + ÿ) 

D’où l’on tire la valeur moyenne sur une période : 
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1 ' 9 

< Pd >= 2 aX o 

L’étude des variations de la valeur moyenne de la puissance < P >=< Pp >=< Pd > en 
fonction de la pulsation d’excitation montre que la valeur maximale de la puissance moyenne 
est obtenue pour fi = co 0 quelle que soit la valeur de S. La valeur maximale de la puissance 
moyenne dissipée ou fournie vaut dans ce cas 

F 2 

< P > = LL 

2a 

La figure ci-dessous représente les variations, en fonction de fi, de la puissance moyenne 
dissipée par les forces de frottements ( ou de la puissance moyenne fournie par la force extérieure 

)• 



Courbe de résonance pour la puissance 


f) Bande passante 

On définit par bande passante, la bande des pulsations autour de fi = co 0 pour lesquelles 
< P >>< P > m ax /2. Les deux pulsations fii et fi 2 , situées de part et d’autre de la pulsation 
ujo et pour lesquelles < P >=< P > max /2, sont appelées pulsations de coupure. La bande 
passante B s’écrit : 


B = fi 2 — fii 

Le calcul de B consiste à rechercher les deux pulsations pour lesquelles < P >=< P > max / 2. 
On obtient l’expression de la bande passante B : 
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B = n 2 -n 1 = 2S 

g) Coefficient de qualité d’un oscillateur 

Le coefficient de qualité d’un oscillateur est défini par le rapport de la pulsation propre co 0 
à la largeur de bande B : 


Q = 


U 0 
~B 


3.3.2 Cas d’une excitation périodique 


Nous avons étudié dans le paragraphe précédent la réponse d’un système vibratoire à une 
excitation sinusoïdale dite excitation harmonique. En pratique, les excitations mécaniques ne 
sont pas toujours parfaitement sinusoïdales ; elles sont souvent périodiques. En considérant le 
cas d’excitations périodiques, nous procéderons à une généralisation du cas harmonique. 

Soit une excitation périodique appliquée à un système amorti à un degré de liberté. L’équa- 
tion différentielle qui régit ce système s’écrit : 

q + 2Sq + culq = A(t ) 


La fonction A(t) étant périodique, de période T, son développement de Fourier s’écrit : 


m 


Cio 

T 


oo 

+ a n cos (nut) + b n sin(ncuf) 

n= 1 


L’équation différentielle s’écrit alors : 


q + 2 5 q + cul C 1 


a o 

~2 


OO 

+ a n cos (nut) + b n sin(ncuf) 

n = 1 


La réponse permanente ( ou stationnaire ) qui s’identifie avec la solution particulière, pour 
t suffisamment élevé, peut alors être calculée pour chacune des composantes de l’excitation : 
o 0 /2, a n cos(ncuf), b n sm(nujt). On obtient alors par superposition : 


q(t) 


Clo 

2^ 0 


n = 1 


a n cos {u n t + ^n) + bn sin (u n t + ipj 



AS 2 cu 


2 

n 


3.4 Impédance mécanique 

3.4.1 Définition 

Considérons un système mécanique soumis à une force sinusoïdale F (t) = F () cos (At). En 
régime permanent, le point d’application de cette force se déplace avec une vitesse v ( t ) = 
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Vo cos (Qt + cf)) . On appelle impédance mécanique d’entrée du système mécanique, le rapport 
des amplitudes complexes de la force F et de la vitesse v 


3.4.2 Impédances mécaniques 

Amortisseur 

Dans le cas d’un amortisseur, la force appliquée est reliée à la vitesse par 

F = av 

On en déduit l’impédance complexe d’un amortisseur 

Z_ a = a 


Masse 

Dans le cas d’une masse, la relation fondamentale de la dynamique s’écrit 

clv 

r = m— 
dt 

On en déduit l’impédance complexe d’une masse 

Z_ m = jrrifl = mQ e 7 2 

Ressort 

Dans le cas d’un ressort de raideur k, la force appliquée / appliquée au ressort s’exprime en 
fonction de l’allongement par 

f = kx 

On en déduit l’impédance complexe d’un ressort 

k .k k _a 7 r 

3.4.3 Puissance 

La valeur moyenne, sur une période, de la puissance fournie est 


< Pf >= 7^F 0 X q cos (0) = ^ Re (Z E ) X { 
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3.4.4 Applications 

Système mécanique résonant 

Soit un système mécanique constitué d’un ressort de raideur k, d’un amortisseur de coef- 
ficient de frottement visqueux a et d’une masse m soumise à une force sinusoïdale F (t) = 
F 0 cos (fit). L’impédance d’entrée de ce système est 


Z E = a+j 



k 

n 


A la résonance fil = cuq = 


, le module de l’impédance est Z E = a. Lorsque la pulsa- 


tion fil — > oo, l’impédance Z_ E — jmfl. 




Module de l’impédance d’entrée Amplitude de la vitesse 

Système antirésonant 

Considérons un circuit constitué par un ressort de raideur k dont une extrémité est reliée 
à une masse m et dont l’autre est soumise à une force sinusoïdale F ( t ). Soit x le déplacement 
de la masse m et soit y le déplacement du point d’application de la force F ( t ). Pour calculer 
l’impédance d’entrée de ce système, nous devons d’abord écrire les équations différentielles du 
mouvement : 


mx = k (x — y) 

F = k (x — y) 

En utilisant la notation complexe, on obtient l’impédance d’entrée : 
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& = f = -ï 


km 


mQ - 


La pulsation d’antirésonance est cuo = \ — • Lorsque kl = cuo, la vitesse Y est nulle tandis 

V m 

que le module de l’impédance est oo. Lorsque la pulsation Q — ’ oc. l’impédance Ze — > 0. 




Module de l’impédance d’entrée Amplitude de la vitesse 
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Chapitre 4 

Oscillations libres des systèmes à deux 
degrés de liberté 

4.1 Introduction 

Les systèmes qui nécessitent deux coordonnées indépendantes pour spécifier leurs positions 
sont appelés systèmes à deux degrés de liberté. 

Exemples 





Figure 1 


Figure 2 


Figure 3 


- Figure 1 Si les masses Toi et 7712 sont astreintes à se déplacer verticalement, 2 coordonnées 
x i et a : 2 sont nécessaires pour spécifier la position de chaque masse à chaque instant. 

- Figure 2 Si la masse M est astreinte à se déplacer dans un plan vertical, deux coordonnées 
sont nécessaires pour spécifier la configuration du système. L’une de ces coordonnées 
peut être le déplacement x qui correspond à la translation verticale de la masse. L’autre 
coordonnée peut être le déplacement angulaire 9 pour tenir compte de la rotation de la 
masse. Ces deux coordonnées sont indépendantes l’une de l’autre. 

- Figure 3 Dans le cas du double pendule, deux coordonnées sont nécessaires pour spécifier 
la position des masses m ] et to 2 . Plusieurs choix sont pourtant possibles, en effet on peut 
choisir (a?i, a? 2 ) ou (yi,y 2 ) ou (0i,0 2 )- 
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Il est possible de spécifier la configuration d’un système à l’aide de plusieurs ensembles de 
coordonnées indépendantes ; un ensemble quelconque de ces coordonnées est appelé coordonnées 
généralisées. Il y a autant d’équations de Lagrange que de degrés de liberté ou de coordonnées 
généralisées. Pour l’étude des systèmes à deux degrés de liberté, il est nécessaire d’écrire deux 
équations différentielles du mouvement que l’on peut obtenir à partir des équations de Lagrange 


f cl 

~dc 

dC 

1 clt 

dqi 

dqi 

1 d 

'dC 

dC 

{ Jt 

_dq 2 

dq 2 


4.2 Systèmes à deux degrés de liberté 

4.2.1 Système masses-ressorts en translation 





x 2 


Considérons le système ci-dessus, constitué de deux masses miet m 2 reliées respectivement 
par deux ressorts de raideur k \ et k 2 à deux bâtis fixes. Les deux masses sont reliées par un 
ressort de raideur K. Ce ressort est appelé ressort de couplage. 

Equations différentielles du mouvement 

Les équations du mouvement pour ce système à deux degrés de liberté peuvent être obte- 
nues à partir des équations de Lagrange pour chaque coordonnée Xi(t) et x 2 (t). Soit T et U 
respectivement l’énergie cinétique et l’énergie potentielle 
T = \ m 1 x\ + \m 2 x \ 

U = \ki x\ + \K (aq - x 2 ) 2 + \k 2 x 2 
U = \ (ki + K) x\ + \ (k 2 + K) xl - Kxix 2 
Le lagrangien L = T — U s’écrit alors 

C = i toi ±1 + m 2 x\ - i (À* + K) x{ — ^ (k 2 + K) x\ + Kx rx 2 
Les équation de Lagrange s’écrivent 
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f d 

'dC 

de 

\ dt 

dii 

dx\ 

1 d 

'dC 

de 

[ dt 

_dx 2 

dx 2 


D’où le système d’équations différentielles du mouvement 


f m\X i + (ki + K) x\ — Kx 2 = 0 
\ m 2 âi 2 + (k 2 + A') x 2 - ADq = 0 

Les termes —Kx 2 et —Kxi qui apparaissent respectivement dans la première et la seconde 
équation sont appelés termes de couplage, et les deux équations différentielles sont dites cou- 
plées. 


Résolution des équations différentielles 

Les deux solutions de ces deux équations différentielles sont des fonctions périodiques et sont 
composées de deux fonctions harmoniques de pulsations différentes et d’amplitudes différentes. 
Supposons que l’une de ces composantes harmoniques s’écrive 

xi(t) = Ai cos(euf + 0) 
x 2 (t) = A 2 cos(c ot + (f)) 

où Ai, A 2 et à sont des constantes et eu l’une des pulsations propres du système. La substi- 
tution de x | et x 2 dans le système d’équations différentielles donne 

f [ki + K — mi û; 2 ] A i — K A 2 = 0 
( — K A\ + [k 2 + A — m 2 u j 2 \ A 2 = 0 

Ce qui constitue un système d’équat-ions linéaires homogènes dont les inconnues sont Ai et 
A 2 . Ce système admet une solution non identiquement nulle seulement si le déterminant A (eu) 
des coefficients de A i et /1 2 est égal à zéro. 

[h + K — 777i CU 2 ] — K 

- K [k 2 + K - 777 2 eu 2 ] 

Le déterminant A (eu) est appelé déterminant caractéristique. L’équation A (eu) = 0 est 
appelée l’équat-ion caractéristique ou équation aux pulsations propres. Elle s’écrit 


A (tu) = 


h + K — rriiOü 2 k 2 + K — m 2 uj 2 — K 2 = 0 


ou encore 


4 2 

eu — eu 


h + K , k 2 + K 


77li 


m 2 


kik 2 + ki K + k 2 K 

777 1 777 2 


= 0 


Cette équation est une équation quadratique en eu qui admet deux solutions réelles positives 
eui et eu 2 appelées les pulsations propres du système 
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Cet exemple montre qu’il y a en général deux pulsations propres dans un système à deux 
degrés de liberté. Chacune des coordonnées, aqet x 2 , possède deux composantes harmoniques 
de pulsations uqet cu 2 


xi = A u cos(cüi t + fa) + A 12 cos(c d 2 t + fa) 
x 2 = A 21 cos(cui t + fa) + A 22 cos(c o 2 t + fa) 

où Au, A 12, A 2 i , A 22 , et fa sont des constantes. Le terme de plus basse fréquence 

correspondant à la pulsation oq est appelé le fondamental. L’autre terme, de pulsation cu 2 , est 
appelé harmonique. 

Les doubles indices sont utilisés pour les amplitudes des différentes composantes harmo- 
niques ; le premier indice se réfère à la coordonnée et le second à la pulsation. Par exemple A 12 
est l’amplitude de aq(t) à la pulsation lu 2 . 

Lorsque A 12 = A 22 = 0, x \ et x 2 correspondent à la première solution particulière sont des 
fonctions sinusoïdales, en phase, de pulsation eu i ; on dit que le système oscille dans le premier 
mode. Dans ce cas 


x\ = Au cos(cuif + fa) 
x 2 = A 2 1 cos(cui t + fa) 


Lorsque A n = A 21 = 0, uqet x 2 correspondent à la seconde solution particulière et sont des 
fonctions sinusoïdales, en opposition de phase, de pulsation cu 2 ; on dit que le système oscille 
dans le second mode. Dans ce cas 


Xi = A 12 cos(c ü 2 t + fa) 
x 2 = A 2 2 cos(uj 2 t + fa) 


Etudions les particularités de ces deux solutions particulières : 
La première solution particulière s’écrit : 


Xi = Au cos(cuif + fa) 

X 2 = A 2 1 cos(cup + fa) 

xi et x 2 doivent vérifier le système d’équations différentielles, ce qui donne 

f [h + K - rriiul] A n - K A 21 = 0 
— A Au + [k 2 + A — 7712 <^i] A 2 i = 0 

Ces deux équations permettent d’obtenir le rapport des amplitudes dans le premier mode 
ou fondamental 

A 2 i k\ + K — mi(jj\ K 

^ 1 Au K k 2 + K - m 2 u:\ 
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- La seconde solution particulière s’écrit : 

xi = A 12 cos(c jJ 2 t + <fi 2 ) 
x 2 = A 22 cos(c o 2 t + <j) 2 ) 

x i et x 2 doivent vérifier le système d’équations différentielles, ce qui donne 

f [kl + K — 777-1 <^1] -4-12 — K A 22 = 0 
| ^ — K Ai 2 + [k 2 + K — m 2 l Ü2\ A 22 = 0 

Ces deux équations permettent d’obtenir le rapport des amplitudes dans le second mode 
ou harmonique 

A 22 ki + K — miix>\ K 

Ai 2 K k 2 + K - m 2 ijj\ 

- La solution générale {x i,x 2 ) est une combinaison linéaire de ces deux solutions parti- 
cuières. aqet rc 2 s’écrivent alors 

xi = An cos (ujit + cj)^ + A u cos (oj 2 t + 4> 2 ) 
x 2 = /.q ffn cos (uit + 0i ) + /i 2 ffi 2 cos (ix 2 t + i fr 2 ) 

où An, A 12 , 4>i et b 2 sont des constantes d’intégration dont les valeurs sont fixées par les 
conditions initiales. 

4.2.2 Cas particulier de deux oscillateurs identiques 

Calcul des constantes d’intégration 

Considérons le cas particulier de deux oscillateurs identiques tels que mi = m 2 = m et 
k | = k 2 = k. Dans ce cas les pulsations propres sont respectivement égales à 



Les rapports d’amplit-udes correspondant à ces pulsations sont respectivement // , = +1 et 
Ih = - 1 - 

Soit x io, x 2 q, xm et x 2 o les valeurs initiales respectives de xi, x 2 , xi et x 2 . Tenant compte de 
ces conditions initiales, on obtient le système d’équat-ions suivant qui permet de déterminer les 
constantes d’intégration An, A i2 , (j) l et (j) 2 

An cos^) + A u cos (<j) 2 ) = aqo 
A n cos(0i) - A u cos(0 2 ) = x 20 
-ui An 8111(0^ - u 2 Ai 2 sin (0 2 ) = ii 0 
-oji An 8111(0^ +lü 2 Ai 2 sin(0 2 ) = x 20 

Les solutions de ce système d’équations sont 
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ou encore 


3Û0 + %20 
2 cos^) 


et 


3Û0 — ^20 
2 cos ((f> 2 ) 


Ê10 + ^20 
2 oji sin(0 1 ) 


et 


2I12 — 


È20 ~ &10 

2 uü 2 sin(0 2 ) 


1. Considérons le cas particulier suivant x w = X 20 = .r'o et x w = x 2 q = 0 ; on obtient dans 
ce cas (f) l = (t >2 = 0 , A 12 = Oet A n = x 0 ; d’où 


X\ = x 0 cos(a>it) 
x 2 = cos(a>it) 


Pour ces conditions initiales particulières, les deux masses oscillent en phase à la même 
pulsation uj \ . On dit que le système oscille dans le premier mode. 



2. Considérons un autre cas particulier pour lequel xio = — X 20 = xq et iq 0 = .i'20 = 0 . On 
obtient dans ce cas 4> 1 = </> 2 = 0, A n = 0 et A 12 = x 0 ; d’où 

X\ = Xq COs(cU2 t) 

X 2 = —X 0 COs(u 2 t) 


On dit que le système oscille dans le second mode car les deux masses oscillent en oppo- 
sition de phase avec le même pulsation ut 2 - 
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3. Considérons enfin le cas particulier suivant x w = x 0 , x 2 o = 0 et x 10 = x 20 
0i = 4*2 = 0; ^ii = -4-12 = xq/2. Les solutions s’écrivent alors sous la forme 


0 ; d’où 


x l(t) = cos ( W 1 t) + ^ cos (co 2 t ) 

X 2 (t) = ^ COS (oui t) — COS (cu 2 t) 

Les solutions ne sont plus des fonctions purement sinusoïdales du temps mais des combi- 
naisons linéaires de deux fonctions sinusoïdales de pulsations respectives eu i et uo 2 . x\ et 
x 2 peuvent s’écrire sous la forme 


La figure suivante 
(c’est-à -dire si K 


X\ ( t ) = Xq cos 

x 2 ( t ) = xq sin 


0J 2 - CO 1 
. 2 
C0 2 — Lüi 


t J CO S 
t ) sin 


oo 2 + oui 
v 2 1 
CO 2 + COi 


représente le résultat obtenu dans le cas où coi est très différent de co 2 
» k). 



Si coi est peu différent de co 2 (c’est-à -dire si K « k). on observe un phénomène de bat- 
tement (voir figure ci-dessous). 
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Coordonnées principales 

Considérons les coordonnées p\ et p 2 obtenues à partir des coordonnées x \ et x 2 par les 
relations 


Pi 

P2 


xi + X2 
2 

X\ — X2 
2 


Tenant compte des expressions de x i et x 2 et des valeurs particulières de // , et p 2 pour 
l’exemple étudié, on obtient 


Pl = — cos (uit) 

P 2 = y cos (u 2 t) 

On remarque que, quelles que soient les conditions initiales, p\ et p 2 sont des fonctions 
purement sinusoïdales du temps de pulsations respectives wiet co 2 . Ces coordonnées particulières 
sont appelées coordonnées principales. On peut vérifier que le système d’équations différentielles 
qui régit le mouvement du système considéré s’écrit sous la forme de deux équations découplées 

Pi + u\ pi = 0 
p 2 + üj% p 2 = 0 

Les relations inverses suivantes 


Xl=Pl+ P2 
X 2 =Pl~ P2 

permettent d’obtenir les coordonnées x \ et x 2 à partir des coordonnées principales p\ et p 2 . 
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4.2.3 Pendules couplés 

Considérons le cas de deux pendules simples identiques couplés par un ressort de raideur K 
et qui effectuent des oscillations de faible amplitude repérées par les angles 0 1 et 9 2 . 



Etablissons tout d’abord les équations différentielles du mouvement dans le cas des oscilla- 
tions de faible amplitude. Il est aisé de montrer que l’énergie cinétique et l’énergie potentielle 
s’écrivent sous les formes quadratiques suivantes 

T = ^ml 2 6 1 + ^ml 2 d 2 

U = \ [Kl 2 + mgl \ 9\ + § [Kl 2 + mgl] 9\ - Kl 2 6 x 9 2 

On remarque la présence du terme de couplage ^ Kl 2 6i9 2 dans l’expression de l’énergie 
potentielle. Comme dans l’exemple précédent, on dit que le couplage est élastique. Si le terme 
de couplage n’existe que dans l’expression de l’énergie cinétique, on dit que le couplage est de 
type inertiel. 

Les équations de Lagrange permettent d’obtenir les équations différentielles du mouvement 

ml 2 'Ù l + [Kl 2 + mgl] 9 1 - Kl 2 9 2 = 0 
—Kl 2 9i + ml 2 'Ù 2 + [Kl 2 + mgl] 9 2 = 0 

En l’absence d’amortissement la solution de ce système d’équations différentielles est de la 
forme 


6i(t) = Ai cos(c ut + 4>) 

9 2 (t) = A 2 cos(ujt + 4>) 

Ces deux expressions doivent satisfaire le système d’équat-ions différentielles, d’où 

[Kl 2 + mgl - ml 2 oj 2 } A 1 - Kl 2 A 2 = 0 
—Kl 2 Ai + [Kl 2 + mgl - ml 2 oj 2 } A 2 = 0 

Ce système d’équat-ions admet des solutions non nulles seulement si eu est solution de l’équa- 
tion aux fréquences 
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[Kl 2 + mgl — ml 2 uj 2 ~\ 2 — K 2 l 4 = 0 
D’où l’on tire l’expression des pulsations propres uqet cu 2 



Les solutions du système d’équations différentielles sont donc 

0 1 = An cos(cuif + (f>i) + Ai 2 cos(cc>2 1 + 0 2 ) 

02 = A 2 1 COS (u>i t + (pi) + A22 cos (u 2 t + (j) 2 ) 

Pour calculer les rapports des amplitudes dans les modes, on suppose que le système oscille 
soit dans le premier mode soit dans le second mode. Dans le premier mode, on obtient le système 

[Kl 2 + mgl — ml 2 üjf\ — Kl 2 ^ = 0 
—Kl 2 + [Kl 2 + mgl — ml 2 cul] /.q = 0 

Dans le second mode, on obtient 

[Kl 2 + mgl — ml 2 uj^\ — Kl 2 g 2 = 0 
—Kl 2 + [Kl 2 + mgl — ml 2 cu|] /i 2 = 0 

Tenant compte des expressions de uj 1 et to 2 on obtient les valeurs du rapport des amplitudes 
dans les modes /i 1 = +1 et /j 2 = —1. Les solutions du système d’équations différentielles 
s’écrivent alors 


0 1 — An cos ( Uit + (f>i) + Ai 2 cos [uj 2 t T <^2) 

0 2 = A n COS (üJit + (fii) - A12 cos (oj 2 t + (fi 2 ) 
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Oscillations forcées des systèmes à 
deux degrés de liberté 

5.1 Equations de Lagrange 

Soit- un système à deux degrés de liberté, soumis à des forces qui dérivent d’un potentiel, à 
des forces de frottement de viscosité et des forces extérieures. Si les coordonnées généralisées 
sont q 1 et q 2 , les équations de Lagrange s’écrivent : 

d_ \dC] _d£ L àV__ 

dt dqi dqi ^ dqi qi 

d_ \dC] _d£ àV__ 

dt _dq 2 _ dq 2 + dq 2 92 

Dans cette expression F qi et F, n sont les forces généralisées conjuguées des coordonnées 
généralisées respectives q\ et q 2 . Elles sont respectivement définies par 

- F qi = 5qi ^ 0 , dans cette expression èlf j représente le travail des forces extérieures pour 

5 ^ 2=0 

une variation Sqi de la coordonnée q\ . lorsque 5q 2 = 0. 

- F q2 = 5(J1=0 , dans cette expression SW 2 représente le travail des forces extérieures pour 

Sq2^=0 

une variation ôq 2 de la coordonnée q 2 , lorsque Sq 2 = 0. 

5.2 Système masses-ressorts-amortisseurs 

Pour étudier les particularités des oscillations forcées des systèmes à deux degrés de liberté, 
étudions le système symétrique suivant soumis à une force horizontale F, appliquée à la première 
masse. 
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5.2.1 Equations différentielles 

Les équations différentielles du mouvement s’écrivent : 

mx i + (k + K) Xi + axi — Kx 2 = F 
—Kx i + mx 2 + (k + A") 0:2 + œc 2 = 0 

5.2.2 Etude du régime permanent sinusoïdal 

Solution permanente 

La solution générale de système d’équations différentielles est égale à la solution de la solu- 
tion du système homogène et d’une solution particulière. La solution de l’équation homogène, 
en raison de l’amortissement, tend vers zéro lorsque le temps augmente. Lorsque le régime 
permanent s’établit, la solution devient égale à la solution permanente et s’écrit alors : 

Xi = Xi cos (fit + q iq) 
x 2 = X 2 cos (Qt + é 2 ) 

Pour calculer les amplitudes X\ et X 2 , ainsi que les phases d, et d 2 . utilisons la méthodes 
des nombres complexes. On peut ainsi écrire : 

Xj. = Re (X, e jüt ) a ; 2 = Re (X 2 e jüt ) F = Re ( F e jQt ) 

Dans ces expressions les amplitudes complexes sont définies par 

Xi = Xi e j ^ X 2 = X 2 e j ^ F = F 0 e j0 

Dans ce cas les équations différentielles se transforment en équations algébriques : 

f (k + K - mO 2 + jaÜ) X x - KX 2 = F 
\ -KX, +(k + K - mO 2 + jaÜ) X 2 = 0 
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Amortissement négligeable 


Considérons d’abord le cas d’un amortissement suffisamment faible pour que l’on puisse 
considérer que a ~ 0. Le système d’équations différentielles s’écrit alors 


f (k + K - mfi 2 ) X x - KX 2 = F 
1 ~KX i + (k + K - mfi 2 ) V 2 = 0 


Les solutions de ce système sont : 


F (Ü 2 a - fi 2 ) 

m (uj\ — fi 2 ) (cu 2 — fi 2 ) 


= 


KF 


m * [üj 


î 


, . k k+2 K 

Les pulsations uj\ = \ — et cu 2 = \ 

V m V m 

précédent. La valeur de la pulsation fi,i est : 


1 

- fi 2 ) (a; 2 - fi 2 ) 

sont les pulsations propres calculées au chapitre 


fiA 


k + K 


m 


Les amplitudes des déplacements Xi et X 2 sont alors données par 


^2 


f |fii-fi 2 | 

m | uj\ — fi 2 | \cu 2 — fi 2 | 
KF 1 

m 2 \ui\ — fi 2 | |oa| — fi 2 | 


Les variations des amplitudes Xi et X 2 sont représentées sur les figures ci-dessous 



Variation de Xi en fonction de fi 


Variation de X 2 en fonction de fi 


On remarque que le phénomène de résonance se produit pour Xi comme pour X 2 lorsque la 
pulsation d’excitation fi est égale à l’une des pulsations propres u> 1 ou cv 2 du système. L’amortis- 
sement étant très faible, les amplitudes à la résonance sont très importantes. Lorsque la pulsa- 
tion fi devient très grande, ces amplitudes tendent vers zéro. Enfin lorsque fi = fi t, l’amplitude 
Xi est égale à zéro ; pour cette raison, la pulsation fi i est appelée pulsation d’antirésonance. 
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5.3 Impédance 

Considérons le système à deux degrés de liberté étudié dans le paragraphe précédent dans 
lequel nous supposons que l’amortissement est nul (a ~ 0 ). En régime stationnaire, on obtient 
pour l’amplit-ude complexe de la vitesse A , : 

x = : tt F 

On en déduit l’impédance d’entrée : 

v F .m (O 2 — uJf) (O 2 — cu|) 

- E = =J n 

Les figures ci-dessous donnent les variation de X et Z E en fonction de O. On note le phé- 
nomène de résonance lorsque la pulsation d’excitation O est égale à l’une des deux pulsations 
propres uq ou üj 2 . A ces pulsations, le module de l’impédance d’entrée est nul. Enfin, lorsque 
O est égale à la pulsation d’antirésonance O i, la vitesse de la première masse est nulle et le 
module de l’impédance d’entrée est infini. Lorsque O — > 00, Z E ~ mil. 



5.4 Application 

Le phénomène d’antirésonance peut être avantageusement utilisé pour supprimer une vibra- 
tion résultant d’une résonance dans un système mécanique. 

Considérons le système à deux degrés de liberté de la figure ci-dessous. 
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Deux degrés de liberté Variation de Xi en fonction de U. 

Les équations différentielles du mouvement s’écrivent 

J mx i + ax i + (k + K) X\ — Kx 2 = F 
\ —Kx i + mâ ;2 + Kx 2 = 0 


En régime permanent sinusoïdal, on obtient 

F 0 


2L i = 


*2 = 


D 2 - — 

“ M 


™ [Sî 4 + n 2 (br + ë) - =] + V - f ) 

ifF„ 1 


Mm [fi- + + f ) - S] + (O 2 - f ) 


Lorsque la pulsation de la force excitatrice est égale à cc.i = \pti-, la niasse m est immobile 

(Xi = 0). 

Si on choisit K et M telles que F = K (c’est-à-dire telles que cu 0 = Qa), la masse m est 

immobile lorsque la pulsation excitatrice O est égale à lü 0 = Dans ces conditions, 

l’ajout de M et K permet d’annuler la vibration de m à cette pulsation. Un tel dispositif 
constitue un "étouffeur” dynamique de vibrations. 
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Chapitre 6 

Généralités sur les phénomènes de 
propagation 


6.1 Propagation à une dimension 


6.1.1 Equation de propagation 

Dans les phénomènes vibratoires traités dans les chapitres précédents, nous nous sommes 
intéressés à des phénomènes ou des grandeurs physiques qui dépendaient d’une seule variable, 
le temps. Nous allons maintenant examiner toute une une série de phénomènes qui sont décrits 
par une fonction qui dépend à la fois du temps t et d’une variable d’espace , x par exemple. 

Ces phénomènes sont régis par une équation aux dérivées partielles, appelée équation d’onde 
ou équation de propagation à une dimension de la forme : 


d 2 s 1 d 2 s 

dx 2 V 2 dt 2 


( 6 . 1 ) 


dans laquelle V est une grandeur physique qui a les dimensions d’une vitesse et sera appelée 
dans la suite vitesse de propagation. 


6.1.2 Solution de l’équation de propagation 

Méthode de D’Alembert 

Pour résoudre l’équation des ondes à une dimension, opérons le changement de variable 
suivant : 


V = (6-2) 

£ = t+^ (6.3) 

Calculons les dérivées partielles par rapport à t et x, en fonction des dérivées partielles par 
rapport à rj et £. 
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Sachant que : 


on obtient 


drj <9£ drj dt; 1 

dt = dt = St qUe âx = ~~Sc = V 


( 6 . 4 ) 


ds ds drj ds dt; ds 

dt drj dt <9£ dt dq 

ds ds drj ds dt; 1 

dx drj dx + dt;dx V 

En tenant compte de ces résultats et sachant que 

d 2 s d 2 s 

drjdi ; dt;drj 


ds 

f âç 

ds ds 
drj dt; 


on obtient : 


d 2 s 

d 2 s 

- 2- 
c 

dF 

drj 2 

d 2 s 

1 

' d 2 s 

dx 2 

V 2 

drj 2 


d 2 s 


d 2 s 

ât2 


- 2 - 


d 2 i 


d 2 s 


( 6 . 5 ) 

(6.6) 

( 6 . 7 ) 

(6.8) 

( 6 . 9 ) 


En remplaçant dans l’équation d’onde et par les expressions ci-dessus, on obtient 
l’équation d’onde exprimée en fonction des dérivées partielles par rapport aux variables r/ et t; : 

d 2 s 


dr/dt; 


= 0 


Cette dernière équation peut s’écrire 


d_ 

drj 


ds 

di 


= 0 


Un intégration par rapport à r/ donne : 


ds 

rr fM 


( 6 . 10 ) 


( 6 . 11 ) 


( 6 . 12 ) 


où / ( rj ) est une fonction qui ne dépend que de rj (et pas de Ç). Enfin une intégration par rapport 
à £ donne : 

s(î|,e) = F(i|)+e($) (6.13) 

où F (rj) , qui ne dépend que de rj , est une primitive de / (rj). La fonction G (£) est une fonction 
qui ne dépend que de £. En revenant aux variables x et t, on obtient la solution générale de 
l’équation des ondes à une dimension : 


s (x, t) = F (t - + G (t + ^ 


( 6 . 14 ) 


Les fonctions F (t — y) et G (t + y) sont des fonctions dont la nature est fixée par les 
conditions aux frontières imposées à la solution s (x, t). 
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Propriétés des solutions particulières F (t — et G (t + |i) 

Propriétés de F (t — |-) On étudie le cas de la solution particulière F (t — y). Pour cela 
on suppose que les conditions aux frontières sont telles que G (t -| |r) est constamment nulle. 
On considère à l’instant- 1\ un point d’abscisse x\ . La valeur de la fonction s en ce point et à cet 
instant est s (x\,ti)- On recherche à un instant t 2 postérieur à t\ {t 2 > C ) la position x 2 d’un 


point pour lequel la valeur de s est la même que la valeur qu’elle avait en x\ 
problème est formulé par l’égalité suivante : 

à l’inst-ant- t\ . Ce 

s(xi,ti) = s(x 2 ,t 2 ) 

(6.15) 

Ce qui se traduit par 

(6.16) 

Cette équation est satisfaite si 


, Xi x 2 

tl ~v- ,2 ~v 

(6.17) 

D’où la valeur de x 2 : 


x 2 = Xi + V (t 2 - ti) 

(6.18) 

Comme t 2 > t±, x 2 est supérieure à x i et ces deux points sont distants de 


x 2 -xi = V (t 2 - h ) 

(6.19) 


F (t — y) correspond à une onde se propageant dans le sens des x croissants (Voir la figure 
ci-dessous). F (t — est appelée onde progressive et cette expression constituera dans la suite 
la définition d’une onde progressive. 



http: / / dj elouah.ifrance.com 








6.1 Propagation à une dimension 


48 


Propriétés de G (t + |i) On étudie le cas de la solution particulière G (t + y). Pour cela 
on suppose que les conditions aux frontières sont telles F (t — y) est constamment nulle. On 
considère à l’instant- ti un point d’abscisse X\ . La valeur de la fonction s en ce point et à cet 
instant est s (xi,ti). On recherche à un instant t 2 postérieur à t\ (t 2 > t\) la position x 2 d’un 
point pour lequel la valeur de s est la même que la valeur en x± à l’instant- t\ . Ce problème est 
formulé par l’égalité suivante : 

s (xytx) = s (x 2 ,t 2 ) ( 6 - 20 ) 

Ce qui se traduit- par 

g (* 1 + f) = g ( < 2 + f) (6 - 21) 

Cett-e équation est satisfaite si 

h + y = h + ^ ( 6 . 22 ) 

D’où la valeur de x 2 : 

x 2 = Xi - V (t 2 - ti) (6.23) 

Comme t 2 > t\, x 2 est inférieure à x\. Ces deux points sont distants de 

xi - x 2 = V (t 2 - ti) (6.24) 

G (t + y) correspond à une onde se propageant dans le sens des x décroissants (Voir la 
figure ci-dessous). G (t + y) correspond à une progressive se propageant dans le sens des x 
décroissant. 



6.1.3 Onde progressive sinusoïdale 

On considère une onde progressive se propageant- dans la direction de l’axe des x, telle que 
le point- d’abscisse x = 0 est- soumis à une vibration sinusoïdale de la forme 

s (x = 0, t) = S 0 cos (eut) (6.25) 
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Le point se trouvant à l’abscisse x > 0 aura la même vibration que celle du point x = 0 
mais avec un retard égal à f? : 

s (. x,t ) = S 0 cos u (t — — j (6.26) 

Cette expression constitue la définition d’une onde progressive sinusoïdale (ou harmonique) ; 
elle peut être écrite sous la forme : 

s(x,t) = So cos [c ot — (f>(x)] (6.27) 

où à (x) = yx représente le déjrhasage lié au temps de propagation y. On dit que (J) (x) 
représente le déphasage dû à la propagation. L’onde progressive sinusoïdale s’écrit sous la forme 
suivante qui permet de mettre en évidence la double périodicité (dans le temps et dans l’espace) : 

M)] 

La quantité T = — est la période temporelle tandis que la quantité À = VT est la longueur 
d’onde qui constitue la période spatiale. On peut vérifier aisément que : 

s(x,t + nT) = s(x,t ) (6.29) 

s (x + n\) = s ( x , t ) (6.30) 


S (x, t ) = S 0 COS 27T 


où n est un nombre entier. 

L’onde progressive s’écrit souvent : 


s (x, t) = So cos [c ot — kx] 

où k = y = ^ est appelé le module du vecteur d’onde qui s’exprime en m 
On utilise très souvent la notation complexe d’une onde progressive sinusoïdale : 


s (x, t ) = So e i{ujt ~ kx) 
s(x,t ) = Se iujt 


(6.31) 


(6.32) 

(6.33) 


où S_ = Sq e~ lkx représente l’amplitude complexe de l’onde progressive sinusoïdale. Le 
module S 0 de S est l’amplitude de l’onde tandis que son argument —kx représente le déphasage 
dû à la propagation. 


6.1.4 Superposition de deux ondes progressives sinusoïdales 

Cas de deux ondes de même fréquence se propageant dans le même sens 

Considérons deux ondes de même fréquence et de même direction de propagation, d’ampli- 
tudes respectives Si et S 2 , et de phases respectives à\ et Ç> 2 - L’onde résultante sera alors : 

s {x, t) = Si e j ( ut - kx+ * i) + S 2 e j{wt - kx+4> = S (6.34) 


http: / / dj elouah.ifrance.com 





6.1 Propagation à une dimension 


50 


ou encore en notation réelle : 


avec 


s (x, t) = S cos ( eut — kx + (j)) 


(6.35) 


S = yj SI + SI + 2 SiSz cos (^ - (j> 2 ) 
et 


(6.36) 


cj) = Arctg 


f S i sin( ( /) 1 ) + S 2 sin(0 2 ) \ 
\Si cos (0j) + S 2 cos (0 2 ) y 


(6.37) 


La superposition de deux ondes harmoniques de même fréquence, et qui se propagent dans la 
même direction, donne une autre onde harmonique progressive de même fréquence, d’amplitude 
S et de phase (j). 


Cas de deux ondes de même fréquence se propageant dans des sens opposés 

Si par contre, on superpose deux ondes harmoniques de même fréquence mais se propageant 
dans des sens opposés , le résultat est tout autre. En effet, dans ce cas : 


s (x, t ) = Sx e ^t- kx +<t>i) + s 2 eii^+kx+M = [S ie j ^e~ jkx + S 2 e j ^e +jkx ] e jujt (6.38) 

et on ne plus écrire l’onde résultante sous la forme d’une onde progressive simple. Un cas 
particulier important se produit quand les deux amplitudes sont identiques. Si on note : 

S 1 =S 2 = S 0 (6.39) 

on a : 


s ( x , t ) 

et donc en notation réelle : 


2^0 cos 






(6.40) 


s (x, t ) = 2S' 0 cos 




cos 



(6.41) 


Ce mode de vibration est très différent d’une onde progressive puisque tous les points x de 
la corde vibrent en phase avec des amplitudes différentes. En particulier, il existe une série de 
points : 


avec 



< fil 02 


27 T 


A 

2 


(6.42) 
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n = 0, ±1, ±2, 

où l’amplitude de vibration est constamment nulle. On dit dans ce cas que l’onde est sta- 
tionnaire et que les points x n sont les nœuds de l’onde. Entre chaque paire de nœuds existe un 
ventre où l’amplitude de vibration est maximum et égale à 2Sq. On note aussi que l’intervalle 
entre deux nœuds est égal à une demi-longueur d’onde À/2. 

6.1.5 Vitesse de phase 

Considérons une onde progressive sinusoidale qui se propage dans le sens des x croissant. 
Un point d’abscisse x possède, à l’instant t, l’élongat-ion : 

s ( x , t ) = S 0 cos (ut — kx) (6.43) 

Entre l’instant t et t + A t l’onde progresse d’une quantité Ax. A l’instant- t + A t, le point 
d’abscisse x + Ax possède la même élongation que celle que possédait le point d’abscisse x à 
l’instant- antérieur t. Ceci se traduit- par l’égalité : 

s ( x , t ) = s (x + Ax, t + A t) (6.44) 

Sq cos (ut — kx) = S 0 cos [u (t + A t) — k (x + Aæ)] (6.45) 

Cette égalité est satisfaite si les phases instantanées sont égales : 

ut — kx = u (t + A t) — k (x + Ax) (6.46) 

Soit- encore 

u At = k Ax (6.47) 

On définit- la vitesse de phase V/ = df qui s’exprime en fonction de u et k apr : 

U = £ (6-48) 

Si la vitesse de phase ne dépend pas de u, le milieu est- dit- non dispersif. Dans le cas contraire 
il est dit dispersif. 

La figure ci-dessous permet d’illustrer la notion de vitesse de phase en considérant- deux 
représentations à des instants différents d’une corde parcourue par une onde . La courbe continue 
représente l’ensemble des points de la corde à l’inst-ant- t. Le point de la corde d’abscisse x est 
représenté par le point blanc, tandis que le point d’abscisse x + Ax est représenté par le point- 
noir. On constate qu’entre les instants t et t + At chacun de ces point- suit- une trajectoire 
rectiligne et- le déplacement du point noir à l’inst-ant- t + At est égal au déplacement du point 
blanc à l’inst-ant- t. En particulier la crête de la corde, correspondant- à une valeur particulière 
de la phase instantanée, semble se déplacer dans le sens de propagation de l’onde avec la vitesse 
de V/, mais la trajectoire de chaque point- matériel est, une trajectoire rectiligne perpendiculaire 
à la direction de propagation. 
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6.1.6 Vitesse de groupe 

La vitesse de phase V n’est pas nécessairement la vitesse que l’on observe lorsqu’on analyse 
un mouvement ondulatoire. En général une onde n’est pas parfaitement sinusoïdale mais a une 
durée limitée et se présente sous la forme d’un train d’onde appelé communément ” puise” ou 
"groupe” qui se propage avec une vitesse Vq appelée vitesse de groupe. Cette onde sous la 
forme d’un puise contient plusieurs fréquences. Si la vitesse de phase est indépendante de la 
fréquence (Milieu non dispersif) alors toutes les fréquences qui constituent le puise se propagent 
à la même vitesse et le puise se propage avec une vitesse de groupe égale à la vitesse de phase. 
Mais si le milieu est dispersif (i.e la vitesse de phase dépend de la fréquence), alors le puise se 
propage avec une vitesse de groupe différente de la vitesse de phase. 

Pour illustrer ce phénomène, considérons une onde constituée de deux ondes de fréquence 
différente et de même amplitude. En x = 0, cette onde s’écrit par exemple sous la forme : 


S (0, t ) = So COS ( Lüit ) + S 0 COS ( U> 2 t ) 

(6.49) 

Cette onde peut s’écrire encore : 


s (0, t ) = 2S'o cos (c o B t) cos (ut) 

(6.50) 

où 

U2 U\ U2 + U i 

u, B = 2 etu,= 2 

(6.51) 


Si üj i est voisine de u 2 - la vibration résultante se présente sous la forme d’une sinusoïde 
de pulsation u dont l’amplitude est modulée par un battement de pulsation u B (Modulation 
d’amplitude). 
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En un point x > 0, l’onde obtenue résulte de la superposition de ces deux ondes qui se sont 
propagées à des vitesses différentes car le milieu de propagation est supposé dispersif : 

s ( x , t) = S 0 cos (cuit — k\x) + Sq cos (u 2 t — k 2 x ) (6.52) 

s ( x , t) peut s’écrire : 


s (x, t ) = 2S' 0 cos (c u B t — ksx) cos (eut — kx ) 
Dans cette expression : 


(6.53) 


kn = 


k2 — k\ h: 2 + ki 

et k = 


2 2 (6 ' 54) 
L’amplitude du battement se propage à une vitesse qui est la vitesse de groupe définie par 
la relation : 

U) B CÜ 2 -CÜ 1 du 

G - k B ~ k 2 - kl ~ dk ( ‘ ^ 

Comme u 2 est peu différente de u\. la vitesse de groupe est définie par : 

V a = - 
G dk 

Tandis que la sinusoïde contenue à l’intérieur du battement se propage à la vitesse de phase : 

u 


V,= 


k 


(6.56) 
)hase : 

(6.57) 



Les flèches verticales noires correspondent au maximum des battements qui se propagent à la 
vitesse de groupe. Les flèches verticales blanches correspondent au maximum des vibrations 

qui se propagent à la vitesse de phase. 
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6.1.7 Onde Vectorielle 


Dans ce qui précède, la quantité s ( x , t) représente une grandeur scalaire, mais certains 
phénomènes décrits par des vecteurs conduisent à des équations similaires : 

d 2 Â 1 d 2 Â 


dx 2 V 2 dt 2 


Le vecteur A, défini dans un milieu trois dimensions, a trois composantes A x . A y . A z et l’ex- 
pression ci-dessus signifie que chacune de ces composantes satisfait individuellement l’équation 
de propagation : 

d 2 A T _ J_ _ n 

dx 2 V 2 dt 2 u 

#-^^ = 0 (6-59) 

d i A z 1_ a 2 A z _ n 

dx 2 V 2 dt 2 u 

Chacune de ces composantes A x , A y . A z se propage en t — y et t + y. 


(6.58) 


6.2 Propagation en trois dimensions 


6.2.1 Equation de propagation 

Dans un système de coordonnées cartésiennes, l’équation de propagation en trois dimensions 
s’écrit sous la forme : 

3 2 s d 2 s d 2 s 1 d 2 s n 

faÂ + ëhj 2 + ]hÂ~V 2 dt 2 ~ ' ' ^ 

On définit le laplacien scalaire de s par l’expression ci-dessous : 


d 2 s d 2 s d 2 s 
S = fa 2 + d^ 2 + fa 2 

et l’équation des ondes s’écrit sous la forme condensée : 

1 d 2 s 


As 


= 0 


(6.61) 


(6.62) 


V 2 dt 2 

s est fonction du temps mais également des cordonnées du point M où la fonction s doit 
être calculée. Si l’on appelle le vecteur position r = OM . la quantité s dépend du temps et du 
vecteur position r ; on écrit s (r,t). 


6.2.2 Onde plane progressive sinusoïdale 

Définition 

L’onde progressive sinusoïdale (ou harmonique), se propageant dans une direction donnée 
par un vecteur unitaire u est définie par : 

s (r, t ) = Sq cos (ujt — k ■ (6.63) 
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où le vecteur k = k ü est appelé le vecteur d’onde. Si les composantes du vecteur k sont k x . 
k y et k z , alors l’onde plane est définie par 

s ( x , y, z, t ) = S 0 cos ( ut — k x x — k y y — k z z) (6.64) 

On peut utiliser la notation complexe pour représenter l’onde plane progressive sinu- 
soïdale qui s’écrit dans ce cas sous la forme : 

s(f,t) = S 0 e< ut -^) (6.65) 

Relation de dispersion 

En remplaçant s(r. t) par son expression dans l’équat-ion de propagation, on obtient la 
relation k = k (eu) pour que l’onde plane définie ci-dessus constitue une solution particulière de 
l’équation d’onde. Cette relation est appelée la relation de dispersion et elle s’écrit : 

k = — ( 6 . 66 ) 

Surface d’onde 

On appelle surface d’onde ou surface équiphase, l’ensemble des points de l’espace pour 
lesquels, au même instant, s (r, t ) a la meme valeur. Recherchons la surface d’onde passant par 
un point M 0 à un instant t ; cette surface est l’ensemble des points M de l’espace pour lesquels 
l’égalité suivante est satisfaite : 


s(r,t ) = s(f 0 ,t) 

(6.67) 

Cette égalité se traduit par : 


S 0 = S 0 

(6.68) 

Cette égalité est satisfaite si 

k ■ r = k ■ fo 

(6.69) 

Tenant compte des propriétés du produit scalaire, on obtient 


k ■ (r — f 0 ) = 0 

(6.70) 


La surface d’onde passant, à l’instant- t, par le point M 0 est l’ensemble des points M satis- 
faisant l’équat-ion ci-dessus. Cette surface est un plan passant par M 0 , et perpendiculaire à la 
direction du vecteur d’onde, donc à la direction de propagation. L’onde est dite plane. 
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Il existe d’autres types d’ondes définis par les surfaces d’onde respectives : par exemple les 
ondes sphériques pour lesquelles les surfaces d’onde sont des sphères ou les ondes cylindriques 
pour lesquelles les surfaces d’onde sont des cylindres. 

Polarisation 

Dans le cas d’une onde plane progressive sinusoïdale représentée par une quantité vectorielle 
A ( f,t ), cette quantité peut avoir différentes orientations par rapport aux surfaces d’ondes : 

1. A est constamment perpendiculaire à la surface d’onde, ou de manière équivalente paral- 
lèle à la direction de propagation : l’onde est dite longitudinale. 

2. A est contenu dans la surface d’onde, ou de manière équivalente perpendiculaire à la di- 
rection de propagation : l’onde est dite transversale. Dans ce cas, l’extrémité du champ 
vectoriel A peut décrire une trajectoire rectiligne : l’onde transversale est dite à polarisa- 
tion rectiligne. Elle peut décrire une trajectoire circulaire (onde transversale à polarisation 
circulaire), ou une trajectoire elliptique (onde transversale à polarisation elliptique). 
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Chapitre 7 

Cordes vibrantes 


7.1 Equation des ondes 

Considérons une corde tendue, rectiligne selon la coordonnée x, et de longueur infinie. Nous 
allons étudier la propagation d’un faible ébranlement le long de la corde. Supposons que cet 
ébranlement se produise suivant l’axe Oy. 

Etudions l’équation du mouvement de cette corde. Nous dénoterons par T la tension à 
laquelle est soumise la corde. On considère en un point d’abscisse x un segment très court de 
cette corde, de longueur Ax. La masse Am du segment est donnée par : 

Am = fiAx 

où /i est la densité linéique de masse de la corde, c’est-à-dire la masse par unité de longueur 
qui s’exprime en kg/m. 



Corde vibrant transversalement 


Dans une situation hors équilibre, le segment n’est plus droit, il présente une courbure. Nous 
considérons des mouvements d’oscillation de la corde de petite amplitude 
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u (. x , t) = u ( x , t) e ÿ 

si bien que nous pouvons faire l’approximation : 


sm 


, du 

= tanWI ' = sï 


sm 


æ+Aæ 


= tan 


x-\-Ax 


du 

dx 


x-\-Ax 


Cette approximation néglige aussi l’allongement du segment, et considère donc la tension T 
comme constante. La force appliquée sur le segment dans la direction y est la résultante de la 
force appliquée au point x , qui est une force appliquée vers le bas et égale en module à 


F (x, t) = T sin 


= T tan 


du 

dx 


et de la force appliquée au point x+Ax qui est vers le haut et égale à 


F (; x + Ax, t) = T sin 


x-\-Ax 


= T tan 


x-\-Ax 


du 

dx 


x-\-Ax 


La force totale dans la direction y est donc : 
R = F {x + Ax, t) — F ( x , t) = T 


du 

du 


dx 

x+Ax dx 

X- 




Nous pouvons appliquer maintenant la loi fondamentale de la dynamique au segment Ax. La 
force dans la direction y doit être égale au produit de la masse Am du segment par l’accélération 
de celui-ci. Donc : 


R = Am 
d 2 u 


IL 


dt 2 


= T 


d 2 u 

~dF 

d 2 u 


dx 2 


Si on définit V = 


L qui a la dimension d’une vitesse, on constate que : 


d 2 u 1 d 2 u 


= 0 


dx 2 V 2 dt 2 

qui est l’équation d’ondes de la corde. V qui est la vitesse de propagation de cette onde. 
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7.2 Ondes progressives harmoniques 

7.2.1 Définition 

Une onde progressive harmonique se propageant selon Ox est définie par : 

u ( x , t) = U 0 cos (ut — kx ) 
ou encore en notation complexe 


u 


(x, t) = U 0 e j{ujt ~ kx) 


où k = y = -y est le module du vecteur d’onde, À étant la longueur d ’onde. 


7.2.2 Force en un point 

On appelle force en un point, la projection selon Oy de la force exercée, en ce point, par la 
partie gauche de la corde sur la partie droite : 


F = ~ T ^r 

ox 

Dans le cas d’une onde progressive sinusoïdale, cette relation devient : 

F(x,t) = jkT U 0 e j{ujt ~ kx) 


La vitesse de particules s’écrit : 

l 

û(x,t) = — =juU 0 

On constate que pour une onde progressive la vitesse de particules ù est en phase avec la 
force F. 


7.2.3 Impédance 

On appelle impédance en un point le rapport de l’amplit-ude complexe de la force à l’am- 
pliude complexe de la vitesse de particule 

Z (*) = J 

Uy 

Dans le cas d’une onde progressive, on obtient : 

Z (x) = = \J )iT 

La quantité \fjFT définit l’impédance caractéristique de la corde 

Z c = s/JÎT = ,jV 
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On obtient une propriété de l’onde progressive plane 

Z (. x ) = Z c Wx 

7.3 Oscillations libres d’une corde de longueur finie 



Corde de longueur L fixée aux extrémités 


Considérons une corde de longueur L fixe aux points x = 0 et x = L. Recherchons une 
solution de l’équation d’onde sous la forme : 


u(x,t) =g(x)f(t) 

En remplaçant dans l’équation de propagation, on obtient : 

1 d 2 g 1 1 d 2 f 

~gd^z = v^Tæ 

Le membre de gauche de cette équation ne dépend que de x. tandis que le membre de gauche 
ne dépend que de t. Ces deux expressions sont donc égales à une constante qui doit être un 
nombre réel négatif que nous posons égal à — k 2 car la solution ne doit pas tendre vers l’infini 
lorsque t tend vers l’infini. Posons lu = kV . On en déduit que : 


d 2 g __ 

d â- 

dt 2 


-h 2 g 


Les solutions de ces deux équations différentielles sont de la forme : 


f = A cos (eut) + B sin (c ut) 
g = C cos (kx) + D sin (kx) 
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La solution de l’équation d’onde peut alors s’écrire sous la forme : 

u (. x , t) = [A cos (c ut) + B sin (eut)] [C cos (kx) + D sin (kx)] 
Tenant compte des conditions aux limites 

u (0, t) = 0 
u ( L , t) = 0 

on obtient 


<7 = 0 

k = nj- o n = 0, 1, 2, 

La solution de l’équat-ion d’onde qui satisfait ces conditions aux limites est donc une somme 
d’une infinités de termes : 


avec 


OO 

u ( x , t ) = [a n cos (cu n t) + b n sin (u n t)] sin (k n x) 

n = 0 


k n = n— et oj n = k n V = n— 

Les u n sont les pulsations propres. Les coefficients a n et b n sont déterminés par les conditions 
initiales du mouvement. Supposons qu’à t = 0 nous imposions à la corde une certaine forme 
initiale u(x, 0) = u 0 (x) et une vitesse initiale 


û ( 0 , t ) = v 0 ( t ) 

Dans ce cas nous aurons les conditions initiales suivantes : 

OO 

u o (t) = a n si n (Kx) 

n = 0 
oo 

(x) = -u n b n sin (. k n x ) 

n = 0 

On doit inverser ces équations pour obtenir les coefficients a n et b n . La méthode de Fourier 
consiste à les multiplier par sin ( k m x) et les intégrer entre 0 et L. Si on utilise les intégrales : 

{ 0 si m n 

L 

— si nn = n 

on obtient 



Cl n 

bn 


- «o (x ) sm 


Un B 


v 0 [X 



sin 



dx 
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7.4 Réflexion et transmission 

7.4.1 Réflexion et transmission entre deux cordes semi-infinies 



incidente réfléchie t 


transmise 



Z 2 - 


Réflexion transmission dans deux cordes semi-infinies 


Soit deux cordes de longueur semi-infinie, reliées en x = 0. Leurs masses linéiques sont 
respectivement /q et /j 2 . Lorsqu’une onde venant de —oc se propage vers x = 0 dans la première 
corde, elle donne naissance au point de jonction, x = 0, à une onde réfléchie et une onde 
transmise. L’écriture de la continuité du déplacement et de la force en x = 0 permet d’obtenir 
le coefficient de réflexion R u et le coefficient de transmission T u définis respectivement par : 


„ _Ur 
““ U, 

T — — 

u ~ Ui 

où Ui , U a et Ut sont les amplitudes des déplacements associés respectivement à l’onde 
incidente, l’onde réfléchie et l’onde transmise. On en déduit : 


R 


U 


Z\ — 2i2 
Z\ + Z2 


T 

-L 11 


2Zi 

Z\ + Z 2 


7.4.2 Réflexion sur une impédance quelconque 
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- fin 


Corde semi-infinie terminée par une impédance Z t 


Soit une corde de longueur semi-infinie, de masse linéique // . tendue horizontalement avec une 
tension T et terminée en x = 0 par une impédance mécanique Zj. Lorsqu’une onde harmonique 
se propage dans la corde de — oo vers x = 0. elle subit une réflexion en ce point. Sachant que le 
déplacement de particules s’écrit : 


u ( x , t ) = Ui e j ^- kx) + U R e j{ -“ t+kx 1 


on en déduit la vitesse de particules et la force en un point d’abscisse x 


ù=^= ju [Ui e ^ ut ~ kx ) + U R e ^ ut+k U] 

F = -Tg = —jkT [Ui e *(<*-**) _ Ur e j(ut+kx)^ 


En x = 0, les conditions aux limites s’écrivent 


Zr — 


F(0,t) 
û(0, t ) 


On en déduit le coefficient de réflexion R u en fonction de l’impédance caractéristique Z c et 
de l’impédance Z T placée à l’extrémité de la corde : 


R 


U 


Ur 

Ui 


Z c — Zr 

Zc + Zr 
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Chapitre 8 

Ondes acoustiques dans les fluides 

8.1 Introduction 


Les ondes acoustiques sont des ondes élastiques qui se propagent dans les fluides (gaz ou 
liquides). Il est donc possible d’obtenir l’équation d’onde qui régit la propagation des ondes 
planes dans un fluide par la même démarche que celle que nous avons utilisée pour établir 
l’équation de propagation des ondes transversales dans une corde. 

Dans la suite, nous utiliserons les symboles suivants pour étudier l’onde acoustique qui se 
propage suivant l’axe des x : 

x : coordonnée à l’équilibre d’une particule du milieu. 

u x : composante suivant l’axe des x du déplacement de particule par rapport à la position 
d’équilibre. 

p 0 : masse volumique du fluide à l’équilibre 
P : pression instantanée en un point quelconque 
P 0 : pression à l’équilibre 

p = P — P 0 : surpression ou pression acoustique 
c : vitesse de propagation de l’onde 

On entend par particule, un élément de volume contenant des millions de molécules de 
telle sorte qu’il puisse être considéré comme continu, mais toutefois suffisamment petit pour 
que les grandeurs acoustiques comme la pression, la masse volumique et la vitesse de particule 
puissent être considérées comme constantes dans cet élément de volume. Dans ce qui suit, nous 
négligerons les effets de la gravitation de telle sorte que P 0 et p 0 sont uniformes dans tout le 
milieu. On suppose d’autre part que le milieu est homogène, isotrope et parfaitement élastique, 
c’est-à-dire non dissipat-if. 


8.2 Equation d’onde 

Considérons le cas d’une onde plane émise dans un fluide par une membrane vibrante plane. 
Lorsque celle-ci est au repos, la pression dans le fluide est uniforme et égale à P 0 . En se déplaçant, 
par exemple dans le sens des x positifs, la membrane comprime la couche de fluide adjacente. 
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Cette situation est instable : le fluide se détend en comprimant à son tout la tranche voisine. 
L’onde progresse ainsi de proche en proche par une succession de compressions et de détentes. 


P(x) . 

1 

1 

1 

1 

1 

1 

1 

1 

1 

.P(x+dx) 

n 

s 

il 


i 

1 

i 

X 


x+dx 



► 

► 



u x (x) 

u x (x+dx) 



•4 


dx 




Propagation acoustique d’une onde acoustique 


Soit une tranche de fluide de petite épaisseur Ax située à l’abscisse x. Lorsque la pertur- 
bation atteint ce point, les forces agissant sur cette tranche ne s’équilibrent plus et elle se met 
en mouvement. Soit u x (x, t ) le déplacement à l’instant t du plan d’abscisse x. Soit F x ( x , t) et 
F x (x + Ax, t ) les forces agissant sur la tranche de fluide respectivement en x et x + Ax. Ces 
forces s’expriment par : 


F x (x,t ) = S P (x,t) 

F x (x + Ax,t) = —SP(x + Ax,t) 

La résultante de ces deux forces est : 

A F x = F x (x, t ) + F x (x + Ax, t ) 

A F x = —S [P (x + Ax, t) — P (x, t)] 

En faisant un développement en série de Taylor au premier ordre de P (x,t), on obtient : 

dP 

P(x + Ax, t) = P (x, t ) + — — Ax 

ox 

D’où : 


dP 

A F x = -S — 
ox 

Comme P = Pq + p, la force résultant s’exprime par : 


AF x = ~S^Ax 
ox 
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Sous l’action de cette force, la tranche de fluide subit une accélération et en écrivant la 
relation fondamentale de la dynamique, on obtient : 


Am 


d 2 u x 

dt 2 


-S d l A, 
dx 


Le fluide étant compressible, le déplacement du plan d’abscisse x + Ax est différent du 
déplacement du plan d’abscisse x et il vaut u x ( x + Ax. t). De nouveau un développement en 
série de Taylor au premier ordre permet d’écrire : 

, \ \ / \ du x 

u x (x + Ax, t) = u x (x , t) + — — Ax 

ox 

Pour prendre en compte la compressibilité du fluide, calculons la dilatation volumique subies 
par la tranche de fluide . Soit Avq = S Ax, le volume à l’équilibre et soit Av, le volume en 
cours de mouvement, avec : 


Av 

Av 

Av 

Av 


S [x + Ax + u x (x + Ax) — x — u x (x)] 
S [Ax + u x (x + Ax) — u x (x)] 
dm 


S 
Av 0 


Ax + 


du a 

dx 


dx 

■Aun 


Ax 


On en déduit la dilatation volumique 


Av — Avq 
Avq 


du 

Rappelons que pour un fluide compressible, la surpression p est reliée à la dilatation volu- 
mique 0 par la relation 

p = —k 6 

où k est le module de compressibilité. On obtient ainsi : 


p = —K 


du x 

dx 


Remarque : On utilise souvent le coefficient de compressibilité x 
Les deux équations 


1 

K 
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d 2 u x dp 

Po ~dP = ~&c 

du x 

P = -k — 
dx 

constituent les deux équations fondamentales de l’acoustique. 

En remplaçant dans la première équation p par son expression tirée de la seconde équation 
on obtient l’équation de propagation : 


d 2 p 1 d 2 p 

dx 2 V 2 dt 2 

où V représente la vitesse de propagation définie par 

V = 

V Po Vpôx 


8.3 Vitesse du son 

Le phénomène de propagation étant un processus adiabatique, la relation liant la pression 
et le volume est 

Pu 7 = cons tan te 

En calculant la différentielle, on obtient : 

u 7 clP + Pv 1 - 1 clv = 0 


Si l’on considère que dP représente la varaiat-ion de pression au voisinage de la pression à 
l’équilibre i y. on obtient : 

vl p + 7P0 Av = 0 


d’où : 


p = -7P0 


An 

^0 


En tenant compte de la définition du module de compressibilité, on obtient : 

1 


k = - = 7P0 
X 


D’où la vitesse du son dans un fluide : 


V = 


h p o 

Po 


Exemple : Dans l’air, dans les conditions normales T = 20°C et P 0 = 10 5 iV ■ m 2 , 7 = 1.4 
et Po = 1.29 kg ■ m~ 3 , on en déduit V ~ 330 m ■ s 
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La valeur de la pression à l’équilibre dépend fortement de la température. Pour une mole 
de gaz parfait, on a : 

P 0 v 0 = RT 


/ 7 RT 


Le produit p 0 v 0 représente la masse molaire M du gaz ; d’où : 


l'y RT 


Dans un gaz parfait, la vitesse de propagation du son est proportionnelle à la racine carrée 
de la température mesurée en °K. 

8.4 Onde progressive sinusoïdale 

8.4.1 Définition 

Une onde acoustique sinusoïdale, s’écrit : 


P (x, t) = po cos U (t - -^) 


On définit le module du vecteur d’onde k par 


k~- 

k ~v 


p ( x , t ) = po cos (c ot — kx ) 

En notation complexe, l’onde progressive sinusoïdale s’écrit 

p ( x , t) = po e j{ujt ~ kx) 

La relation liant la pression acoustique et la compressibilité, à savoir 


P = — K 
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permet d’écrire 


u(x, t ) = — p(x,t ) dx 


K 


U 


(x,t) = -- f p 0 e j{ujt ~ kx) dx 


u(x,t ) = 

jkn 

u(x,t) = -J^- e ^ ut - kx) 
jujp 0 V 

La dérivation de cette dernière expression par rapport au temps permet d’obtenir la vitesse 
de particules : 

«(*.*) = Ï = po 


dt p 0 V J 

On constate que pour une onde progressive la vitesse de particules est en phase avec la 
pression acoustique. 


8.4.2 Impédance acoustique 

On appelle impédance acoustique en un point le rapport de l’amplitude complexe de la 
pression à l’ampliude complexe de la vitesse de particule 


Dans le cas d’une onde progressive, on obtient : 

Z (x) = p 0 V 

Le produit p 0 V définit l’impédance acoustique caractéristique du fluide 

Z c = p 0 V 

On obtient une propriété de l’onde plane progressive : 

Z (æ) = Z c Wx 

8.4.3 Energie acoustique 

Densité d’énergie cinétique 

Soit un petit élément de volume vq dont le dépélacement est u (x. t ) et dont la vitesse est 
û(x,t); il possède une énergie cinétique 

P 1 -2 

E c = l^Po v O u 
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On définit l’énergie cinétique par unité de volume ou densité d’énergie cinétique 

„ E c 1 , 

£c = — = -tzPqU 

vq 2 

Densité d’énergie potentielle 

Soit un petit élément de volume vq. Sous l’action de la surpression p, cet élément se comprime 
ou se dilate en raison de la compressibilité du fluide. L’énergie potentielle emmagasinée est égale 
au travail fourni par la pression pour comprimer ou dilater le volume vo : 

E p = [ —pdv 


Sachant que 
on en déduit que : 


p = —K 


v -vp 

Vo 


dv = Vq dp 

K 

d’où : 


E p = — pdp 

K J o 

P V ° 2 

E 'P = ttP 

2 K 

On en déduit la densité d’énergie potentielle : 

g Ep 1 2 

o P 

Vq ÀK 

Densité d’énergie 


La densité d’énergie est égale à la somme de la densité d’énergie cinétique et de la densité 
d’énergie potentielle 


S — S c + Ep 

1 9 1 9 

Dans le cas particulier d’une onde plane progressive sinusoïdale, 

1 


£ = E = 
c p 2p 0 V 2 


Pq cos w (c vt — kx) 
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Tl 


et 


£ = 


1 




Po cos 2 (ut — kx) 


On définit la valeur moyenne temporelle de la densité d’énergie comme 

r t+T 


1 f 

(£} = t [ 


£ dt 


2vr 

où T = — . 
u 

On définit également la moyenne spatiale de la densité d’énergie : 


(£) = 


rx+X 


£ dx 


Dans le cas d’une onde progressive, ces deux valeurs sont égales et valent : 


{£) = 


pI 


2p 0 ^ 2 


Intensité 

On appelle intensité de l’onde acoustique la puissance qui traverse, par unité de temps une 
surface unité perpendiculaire à la direction de propagation. 

Pour calculer l’intensité de l’onde calculons l’énergie qui traverse pendant un intervalle de 
temps une surface S perpendiculaire à la direction de propagation. 


t t+dt 



r 

' 

n 


* 


i 

i 



y k 



i 

i 




> 

s 






U 

i 

L 

il 




◄ ► 

V dt 

Flux de puissance 


Cette énergie dE est égale à l’énergie contenue dans un volume S V' dt et elle égale à 

dE = £ S V dt 
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D’où la puissance V traversant cette surface 


V = = £ SV 

clt 

On en déduit l’expression de l’intensité de l’onde acoustique 

'<*> - § p 
I(t) = £ V 


Ht) = 


Po _ _ _2 


PoV 


cos 2 (c ut — kx) 


On appelle intensité de l’onde acoustique la valeur moyenne 


I 

I 

I 


1 

T J t 

Po 

2p 0 V 

Po 


rt+T 


I ( t ) dt 


2Z r 



Niveau sonore 

On définit le niveau sonore en décibels par 

N dB = lOlog 

I 0 est une intensité de référence correspondant à I 0 = 10 12 W ■ m 2 . 

Exemple : Pour une fréquence / = 1kHz, le seuil d’audition est égal à 0 dB et le seuil de 
douleur est égal à 130 dB. Pour calculer l’intensité, l’amplitude de pression p 0 la vitesse de 
particules et le déplacement de particules, dans le cas où V = 330 m/s et Z c = 411 rayleighs, 
on peut utiliser les relations suivantes : 

I = J 0 10 017Vds 

Po = V2 Z C I 
Po 

Z c 
û 


u = 


u = 


2-nf 


Pour chacun de ces deux cas, on obtient : 


^dB 

I iW/rri 2 ) 

Po {P a) 

û (m/s) 

u (m) 

0 dB 

10" 12 

2. 9 x 10~ 5 

7 x 10~ 8 

1. 1 x 10“ n 

130 dB 

10.0 

91 

0. 22 

3. 5 x 10" 5 


http:/ / dj elouah.ifrance.com 





8.5 Reflexion- Transmission 


73 


8.5 Reflexion- Transmission 


X 


•> 


Onde incidente 



Onde réfléchie 

Pl’ M P2’ ^2 

Réflexion à un interface fluide-fluide 


-► 

X 


Soit deux milieux fluides semi-infinis séparés par un surface plane. Choisissons un repère 
orthonormé de telle sorte que le plan yOz coïncide avec la surface de séparation. Lorsque une 
onde acoustique provenant de — oo, se propageant dans le premier dans la direction de l’axe des 
x arrive à la surface de séparation, elle donne naissance à deux ondes 

- une onde réfléchie qui se propage dans le premier milieu dans le sens des x décroissant. 

- une onde transmise qui se propage dans le second milieu dans le sens des x croissant. 

L’onde résultante dans le premier milieu (x < 0) est caractérisée par : 


Pi(x,t) = Pi(x,t)+p R (x,t) 

Pi(x,t) = p i é>^ t - klX ^ + P R é i( - ut+ kix) 

Ù! ( x , t) = [Pi e j{ujt ~ kix) - P R e j( -“ t+kix) ] 

Zi 

Dans le deuxième milieu, on a 

P2 ( x , t ) = pt (x, t ) 

p 2 (x,t) = P T e j ^- k2x) 

u 2 (x,t) = ±-P T e?( ut - kax '> 

Les relations de continuité à l’interface s’écrivent 

f pi(0,t) =p 2 (0,t) 

\ Ùi (0, t) = ù 2 (0 ,t) 

On en déduit 
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ou encore 


1 

~Zx 


Pi + Pr = Pt 
(P - Pr) = P Pt 

Zj2 


1 


Pi 


R 


l-^ = 

Pi 


P T 

\Pi 

Z\ Pt 

Z~2PÏ 


On définit 

- le coefficient de reflexion pour la pression 


Rp 


Pr 

Pi 


- le coefficient de transmission pour la pression 


Tp = 


Pt 

~Pi 


Les deux relations de continuité s’écrivent alors 


1 + Rp — Tp 

i - Rr = §?> 


On en déduit les coefficients de reflexion et de transmission 


Z 2 — Zi 
Z 2 + Z\ 

2Z 2 

Z 2 + Z\ 

En tenant compte des relations p, = Z\ ù t . pp = —Z] v /,> et pp = Z 2 Ùt, on peut calculer les 
coefficients de reflexion etde transmission pour la vitesse de particules et pour le déplacement 
de particules : 


R P = 
T P = 


Rp 

Tù 


Ru = 


Tu = 


Z\ — Z 2 

Z\ + Z 2 
2Z X 

Z\ + Z 2 


En tenant compte des relations I t = Pf /‘2Z\ : Ir = et I T = P£/2Z 2 . on peut calculer 

les coefficients de reflexion et de transmission pour l’intensité acoustique : 
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olt 


ir = \ z 1 -z 2 y 

1% _Z\ + Z 2 

Ij' AZ\Z 2 
li \Z\ + Z 2 
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Chapitre 9 


Propagation d’une onde électrique 
dans une ligne coaxiale 

9.1 Introduction 

Lors de l’étude des circuits électriques, en régime sinusoïdal forcé, on considère habituelle- 
ment que dans une branche quelconque du circuit, le courant électrique est le même en tout 
point et ne dépend que du temps. Cette hypothèse n’est valable que si la longueur d’onde 
associée au courant sinusoïdal est grande devant les dimensions du circuit ; par exemple pour 
une fréquence / = 50Hz, la longueur d’onde correspondante est A = c/ f ~ 6000/cm, où c est la 
vitesse de propagation des ondes électromagnétiques dans le vide (c = 3 x 10 8 m ■ s -1 ). Cette 
hypothèse revient de façon équivalente à négliger le temps de propagation de l’onde électroma- 
gnétique (propagation de vitesse infinie). 

Si À est du même ordre de grandeur des dimensions du circuit, cette approximation n’est 
plus valable ; par exemple si / ~ GHz , À ~ 30cm qui est alors comparable aux dimensions des 
circuits ou des fils couramment utilisés dans les réseaux électriques à l’échelle du laboratoire. 
Dans ce cas, il faut tenir compte de la vitesse finie des ondes électromagnétiques. 

9.2 Equation de propagation 

Un câble possède une capacité par unité de longueur Ci et une autoinductance par unité 
de longueur Li . Nous supposons que ce câble, que nous appelons ligne, est rectiligne et que la 
position de chaque point est repérée par son abscisse x. Soit un élément de câble situé en x et 
de faible longueur dx. Il possède une capacité Ci dx et une inductance L/ dx. Cet élément de 
câble est représenté, à un instant t. par le schéma ci-dessous. 
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4°— ► 

i(x) 

v(x) 

— ■) 

L /dx 

C/dx — 

— ► — O'u 
i(x+dx) 

r 

.U. .(J ^ 

A 




dx 

f 


v(x+dx) 


x 


H ► 

x+dx 


cellule élémentaire 


1. L’écriture des lois de Kirchhoff permet d’établir les relations suivantes : 

dv (. x , t) di ( x , t) 

dx 1 dt 

di (x, t ) dv (x, t ) 

dx 1 dt 

A partir de ces deux relations, on peut établir les équations de propagation de la tension 
v et du courant i : 


Dans ces deux équations, V 
dans le câble . 


d 2 v 1 d 2 v 

dx 2 V 2 dt 2 

d 2 i 1 d 2 i _ 

dx 2 V 2 dt 2 

est la vitesse de propagation de l’onde électrique 


9.3 Solution de l’équation de propagation 


On peut vérifier que les fonctions / {t — y) et- g (t — y) sont- deux solutions de l’équat-ion 
de propagation. La solution de l’équat-ion de propagation s’écrira donc généralement sous la 
forme : 



La fonction / (t — y) correspond à une onde progressive se propageant dans le sens des x 
croissant- tandis que g (t + y) correspond à une onde se propageant- dans le sens des x décrois- 
sant-. 
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9.4 Onde Progressive sinusoïdale 

9.4.1 Définition 

Lorsque g (t + =0, la tension électrique en un point d’abscisse x est donnée par : 

v(x,t) = f(t- 

Dans ce cas, on dit qu’un onde progressive se propage dans le câble. 

Lorsque la fonction f (t — y) est une fonction sinuoïdale, l’onde progressive est dite sinu- 
soïdale ou harmonique : 

V (. X , t ) = V 0 COS U) (t — — j 

On peut encore écrire 

v ( x , t) = vo cos [ut — kx] 

où k = y (Relation de dispersion). 

En utilisant la notation complexe : 

v (x,t) = Re [v 0 e j[ujt ~ kx] ] = Re [v e ju)t ] 

où l’amplitude complexe est définie par : 

v = v 0 e~ jkx 

De même , on peut obtenir l’expression de l’intensité en fonction du temps t et de la position 
x : 

i (x, t) = i 0 cos (ut — kx) 
i (x, t ) = Re [i e jujt ] 

où l’amplit-ude complexe du courant est définie par la relation : 

i = i Q e ~ jkx 


avec io = -%= 

/ Li 


9.4.2 Impédance en un point 

L’impédance en un point est définie par le rapport : 

v 
i 

Dans le cas d’une onde progressive, l’impédance est un nombre réel égal à 

Z c = 

et ne dépend pas de x. Cette impédance est appelée impédance caractéristique ou impédance 
itérative du câble. 



Z{x) 
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Chapitre 10 

Eléments d’analyse vectorielle 


10.1 Champ scalaire - Champ vectoriel 

Soit un trièdre orthonormé (e x ,e y ,e z ) et M un point de l’espace, de coordonnées (x, y. z) : 

> _ 

OM = x e x + y e y + z e z 

La fonction / (M) est dite fonction scalaire de point ou champ scalaire si : 

f ( M ) = f (®, y, z) 

Le vecteur v {M) est dit fonction vectorielle de point ou champ vectoriel si : 

v (M) = V x ( x , y, z ) e x + V y (x, y, z) e y + V z (x, y, z) e z 


10.2 Gradient d’un champ scalaire 


Le gradient ( noté grcicl ) est défini à partir d’une fonction scalaire de point et a pour 
composantes suivant e x ,e y , et e z les dérivées partielles de / (M) par rapport à x, y et z respec- 
tivement : 

— f m - ,df -> ,df _ 

g r ad( f ) = -e x + ^e,+ Ti e 1 


10.3 Divergence d’un champ vectoriel 

La divergence (notée div ) n’est définie qu’à partir d’une fonction vectorielle v (M) de point 
et donne une fonction scalaire de point définie, en coordonnées cartésiennes par : 


div (v) = 


dv 7 


dv y dv z 


dx dy dz 
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10.4 Rotationnel d’un champ vectoriel 


Le rotationnel noté ( rot ) d’un champ vectoriel donne une fonction vectorielle de point 
définie en coordonnées cartésiennes par : 


rot (v) = 


dv z dv y 
dy dz 


+ 


dv T dv~ 


dz dx 


e y + 


dv y _ dvx_ 
dx dy 


10.5 Laplacien scalaire 

Le laplacien scalaire d’une fonction scalaire de point ( noté lap ou A) est par définition un 
champ scalaire défini par : 

lap (/) = A / = cliv grad (/) 

Dans un système de coordonnées cartésiennes, il s’écrit : 

d 2 f d 2 f d 2 f 


lap (/) = A f = 


dx 2 dy 2 dz 2 


10.6 Laplacien vectoriel 


Le laplacien vectoriel (noté lap ou Â) d’un champ vectoriel v est un champ vectoriel défini 
par : 

lap ( v ) = Â (v) = grad [div (u)] — rot rot ( v ) 

Dans le cas d’un système de coordonnées cartésiennes, le laplacien vectoriel a pour compo- 
santes : 


Av r = 


lap ( v ) < 


d 2 v x 

dx 2 
d 2 v 

AVy = ~dâ 

A d 2y Z 

Av ’ = ^ 


d 2 v x 

dy 2 

d\_ 

dy 2 

d\ 

dy 2 


+ 


d 2 v x 

dz 2 

d\_ 

dz 2 

d 2 v z 

dz 2 


10.7 Opérateur nabla 


Pour écrire de manière plus compacte les opérateurs vectoriels précédemment définis, on 
introduit un vecteur symbolique appelé opérateur nabla et défini par : 


d 


d 


d 


^ y x o "b y y o A &z o 

dx dy dz 


Les opérateurs vectoriels s’écrivent parfois à l’aide de l’opérateur nabla sous les formes 
respectives suivantes : 
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le gradient d’un champ scalaire / est noté 


df * df 


df 


grad (/) = V/ = - ë x + - ë y + - e z 


la divergence d’un champ vectoriel est notée 


div (v) = V • v = 


dv x dv v dv 7 

" + - 7 ^ + 


dx dy dz 


le rotationnel d’un champ vectoriel est noté 
rot (v) = V x v = 
le laplacien scalaire d’un champ scalaire est noté 


\dv z dvy 

P, — |— 

dv x dv z 


dv y dv x ~ 

o 

a. 

CO 

— i 

'-'X 1 

dz dx 

1 

dx dy 


lap (/) = A / = div \grad (/) = V • V/ = V 2 (/) 


V 2 se lit ” del de”. 


le laplacien vectoriel d’un champ vectoriel est noté 
lap (■ v ) = Â (v) = grad [ div (u)] — rot \ rot (v) 


= VV ■ (v) - V x 


V x v\ 


10.8 Théorème de Stokes-Théorème de Gauss 


10.8.1 Circulation d’un champ vectoriel 


On définit la circulation d’un vecteur v le 



AB 


long d’un contour (C), par l’intégrale curviligne : 

B 



La circulation de long d’un contour fermé est notée : 
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10.8.2 Flux d’un champ vectoriel 

On définit le flux d’un vecteur v à travers une surface (S) par l’intégrale double 


<t>/{S) {v) = JJ v-nds 

( 5 ) 



Lorsque la surface (S) est fermée, le vecteur unitaire n est dirigé de l’intérieur vers l’extérieur. 


10.8.3 Théorème de Stockes 

La circulation d’un vecteur le long d’un contour fermé (C) limitant une surface (S) est égal 
au flux de son rotationnel à travers cette surface. 


c (v) = 4>/{S) («)) 


' v ■ dl = JJ rot (v) ■ n dS 
(S) 

le vecteur unitaire n est orienté selon la convention du tire-bouchon de Maxwell. 


10.8.4 Théorème de Gauss- Ostrogradski (ou théorème de la diver- 
gence) 

Le flux d’un champ vectoriel à travers une surface fermée (,S')est égal à l’intégrale de sa 
divergence dans le volume (r)limité par la surface fermée (S) : 


v ■ n dS = 



div (v) dr 


( S fermée) 


0 ") 
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Chapitre 11 

Les équations de Maxwell dans le vide 


11.1 Le champ électromagnétique 

Propriétés du champ électrostatique 

Dans les états stationnaires, le champ électrique est appelé champ électrostatique. Le champ 
électrostatique E créé par une distribution de charges de densité p située dans le vide, est à 
circulation conservative, c’est-à-dire qu’il satisfait les relations intégrale et locale : 

j)Ê ■ cil = 0 

(D 

rot (É) = 0 

où (r) est un contour fermé quelconque orienté. 

Il satisfait également satisfait les relations intégrale et locale : 

/ / e ' iS = /// '' 

(S fermée) (r) 

dw ( Ê ) = £ 

où (S) est une surface fermée quelconque orientée vers l’extérieur et (r) est le volume 
intérieur à (5'). e 0 = 8.854187817 x 1CT 12 F m _1 est la permittivité du vide 

11.1.1 Champ électromoteur et vecteur densité de courant 

Le champ électromoteur 

Lorsqu’un courant électrique circule dans un conducteur, cela implique l’existence d’une 
force motionnelle f m agissant sur les porteurs de charge q et l’on définit le champ électromoteur 
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E rn par la relation 




Q 


La circulation de ce champ le long d’un contour fermé orienté (T) n’est pas conservative 
c’est-à-dire qu’elle est différente de zéro. Par définition cette circulation est appelée la force 
électromotrice e relative au contour considéré : 



(T) 


On peut bien entendu définir la f.é.m relative à un tronçon AB orienté, non fermé : 



Le vecteur densité de courant 

Un courant électrique correspond à des charges électriques mobiles. On appelle vecteur 
densité de courant j , le vecteur tangent à la ligne de courant, et défini par 

j = pv 

où p est la densité volumique de charges mobiles et v la vitesse d’entraînement de ces charges 
mobiles. Le module de ce vecteur représente la charge qui traverse par unité de temps, l’unité 
de surface perpendiculaire à la direction de déplacement des charges mobiles ; il s’exprime en 
A.m~ 2 . 

En régime stationnaire, c’est-à-dire lorsque le vecteur densité de courant j est indépendant du 
temps, le flux de j est conservatif ce qui se traduit par les relations intégrale et locale : 

j ■ dS = 0 

(S fermée) 

div (j) = 0 



11.1.2 Le champ magnétique 

Le champ magnétique B créé par une distribution de courants de densité j est à flux 
conservatif, c’est-à-dire qu’il satisfait les relations intégrale et locale : 



div 



0 

0 


http: / / dj elouah.ifrance.com 





11.2 Le régime variable 


85 


où (S) est une surface fermée quelconque. 

Le champ magnétique B satisfait les relations intégrale et locale : 



(U (S) 

rot (b') = Mol 


où (r) est un contour fermé et orienté quelconque et (S) une surface quelconque s’appuyant 
sur (r) et orientée à partir de (r) par la règle dite du "tire-bouchon de Maxwell" ou "du 
bonhomme d’Ampère". Mo = 1-2566370614 x 10 G NA 2 est la perméabilité magnétique du 
vide. 

11.2 Le régime variable 

11.2.1 Le phénomène de propagation 

Considérons un ensemble constitué par des circuits parcourus par des courants et par des 
distributions de charge variant en fonction du temps ; cet ensemble pouvant être au repos ou 
en mouvement. Au voisinage de ces distributions régnent un champ électrique et un champ 
magnétique. Contrairement au cas stationnaire, ces champs ne sont pas synchrones avec les 
sources, c’est-à-dire qu’à un instant t donné, ces champs dépendent des valeurs des sources à 
l’instant t — 9 ; nous exprimons ce fait en disant qu’il y a propagation à vitesse finie des champs 
à partir des sources qui leur donnent naissance et le retard 6 est d’autant plus grand que le 
point où l’on désire connaître les champs est éloigné des sources. 

11.2.2 Le phénomène d’induction 

Un circuit filiforme au repos et parcouru par un courant invariable n’entraîne l’apparition 
d’aucune f.é.m ou d’aucun courant dans un autre circuit filiforme au repos. Il n’en est pas de 
même si le courant varie ou si les circuits en présence se déplacent l’un par rapport à l’autre : la 
f.é.m ou le courant qui apparaissent sont dûs au phénomène d’induction. Ce phénomène entraîne 
l’apparition d’un champ électrique supplémentaire (appelé champ induit) ; ce qui conduit à 
modifier la propriété fondamentale du champ électrique. 

11.2.3 Le phénomène de capacité 

Un circuit comprenant un condensateur alimenté par une source de tension variable en 
fonction du temps, est parcouru par un courant variable bien que la continuité électrique soit 
interrompue par l’espace entre les armatures du condensateur. Dans ce cas l’intensité du courant 
n’est plus conservée tout au long du circuit puisqu’elle est nulle dans l’espace entre les armatures. 
Il n’est donc plus possible d’appliquer le théorème d’Ampère. Pour conserver la validité de ce 
dernier, nous serons amenés à introduire un courant "fictif" appelé courant de déplacement. 


http: / / dj elouah.ifrance.com 





11.3 L’induction électromagnétique 


86 


11.3 L’induction électromagnétique 

11.3.1 Loi de Faraday-Lenz 

On peut induire une f.é.m dans un circuit filiforme ( C ) fermé en faisant varier le flux 
magnétique à travers le circuit : c’est- le phénomène d’induction électromagnétique. 

Pendant un temps dt. la variation du flux magnétique total à travers une surface quelconque 
s’appuyant sur le circuit ( C ) est dé ; la f.é.m induite e s’exprime à l’aide de la loi de Faraday : 

d(f> 

dt 

Cette loi, établie expérimentalement pour des variations relativement lentes du flux magné- 
tique en fonction du temps, est valable pour tout régime variable et elle sert de base à l’étude 
de l’électromagnétisme classique. 

Puisque une f.é.m apparaît dans le circuit (C) et y fait circuler un courant ceci implique 
l’exist-ence d’un champ électromoteur agissant sur les porteurs de charge du circuit (C). Ce 
champ est appelé champ électrique induit. 

11.3.2 Equation de Maxwell-Faraday 

Considérons un circuit ( C ) au repos soumis à un champ variable. Un champ électrique 
va prendre naissance dans tout l’espace où existe un champ magnétique variable. Le champ 
électrique induit joue un rôle de champ électromot-eur et la f.é.m apparaissant dans tout le 
circuit (C) peut s’écrire : 


e — lÉ,.dï=- i 4 = -é 
J dt dt 

(C) 


B -dS 


(S) 


où (S) est une surface orientée s’appuyant sur le contour orienté (C). 
En permutant les opérateurs d’intégration et de dérivation on a : 


f 


(J) Ei ■ dl = 

(C) 


d(B-dS 


dt 


as ,s 

~M' ds 


(S) 


car le circuit (C) étant immobile, la surface (S) l’est aussi et dS est indépendant du temps. En 
appliquant le théorème de Stokes nous pouvons écrire : 

Ei • dl = J Jr°t (4) -dS = -Jj^-dS 

(C) (5) (5) 

Cette égalité étant satisfaite quel que soit (S) s’appuyant sur (C), il en résulte : 

dB 

~dt 


rot ( Ei) = 


http: / / dj elouah.ifrance.com 





11.4 Le théorème d’Ampère 


87 


Remarquons que s’il existe en plus du champ électrique induit un champ électrostatique Es- 
Le champ total E = Ei + Eg satisfait encore la relation ci-dessus car 


rot (j^s) = 0 


et on a la relation de Maxwell-Faraday 

rot (Ê) 


dB 

~dt 


En définitive, nous devons retenir de l’étude du phénomène d’induction électromagnétique 
le résultat fondamental suivant : en chaque point de l’espace où existe un champ magnétique 
variable nous devons associer un champ électrique induit variable à circulation non conservative 
c’est-à-dire ne dérivant pas d’un potentiel. L’ensemble de ces deux champs constitue le champ 
électromagnétique. 


11.4 Le théorème d’Ampère 

11.4.1 Equation de continuité 

Si on considère une surface fermée (S) entourant un volume (r) . Si p est la charge volumique 
et q la charge totale du volume (r) à l’instant- t ; on a 


Q = 



p ch 


(t) 


Pendant- l’intervalle de temps dt. la variation de la charge totale est dq et- on a : 

/ \ 

dq d fff 

T t = jt Jlr dT 

\ M / 


mais 


j-dS représente la charge totale sortant de la surface ( S ) par unité de temps, 


( S fermée) 

donc par suite de la conservation de la charge : 


et 


d_ 

dt 


dq 

dt 


( 


J • dS 



(S fermée) 

\ 


p dr 

V ~ « / 


J • dS 


(S) 
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En changeant l’ordre des opérations par rapport à l’espace et par rapport au temps, on a : 

'dp 



dt 


clr = — 


J • dS 


(r) (S fermée ) 

or, d’après le théorème d’Ostrogradsky : 


j • dS = 


( S fermée) 



div (j\ dr 


donc 



(t) 


% + dtv ( j 


dr = 0 


Cette relation doit être vérifiée quel que soit le volume (r), il faut donc que l’on ait : 

iw ( 7 ) + 1 = 0 

Cette équation dite de continuité traduit la conservation de la charge électrique et montre 
que le flux du vecteur densité de courant n’est plus conservatif comme dans le cas des états 
stationnaires. 

Ce qui se traduit localement par 

div (e) = £- 

v y e 0 

p représente la densité volumique de charges . L’équation de continuité peut alors s’écrire : 


div ( j ) + ^ = div j + ^ (div j = div | j + I = 0 


On voit que la relation fondamentale de conservation du flux de la densité de courant sera 
conservée si on l’applique à une densité de courant total jr égale à la somme de la densité de 
courant vrai j (appelé courant de conduction) et d’une densité de courant fictif appelé courant 
de déplacement et défini par : 


Jd — 


d fe 0 Ê) 


dt 


d’où : 


JT=J+JD=J 


et div (j T ) =0 
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11.4.2 Le théorème d’Ampère 


Le théorème d’Ampère peut être généralisé à condition de l’appliquer au courant total.La 
relation de Maxwell- Ampère qui en est la traduction s’écrit : 

— > / x OÊ 

rot [BJ= /i 0 j +li 0 £ 0 — 

La relation intégrale du théorème d’Ampère généralisé est : 


/ 



ho 




d 



11.5 En résumé 


En dehors des discontinuités, les équations de Maxwell s’écrivent : 
- Théorème de Gauss pour E 
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Chapitre 12 

Propagation des ondes 
électromagnétiques dans le vide 


12.1 Equations de propagation pour E et B 

En absence de charges et de courant, les équations de Maxwell s’écrivent : 


div ( E ) = 0 


( 12 . 1 ) 


div l B) = 0 


rot [E] = — - 


— * ( 

rot [B J = g 0 e 0 — 

Pour établir l’équation relative au champ électrique E , il éliminer le champ magnétique B 
. Pour cela, calculons le rotationnel de chacun des membres de le loi de Fraday : 


rot f rot (£'))= rot — 


(12.2) 


en permutant l’ordre des dérivations, on obtient : 


Sachant que 


— * ( — * E\ <9 ( — » d f dÊ 

rot [rot EJ = — — [rot B J = 1 ^ 0 £ 0 — 


rot rot (É)= grad div ( E) — A Ë 


( 12 . 3 ) 


( 12 . 4 ) 


où A est le laplacien vectoriel, on obtient l’équation aux dérivées partielles suivante : 
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d 

A E - gr ad ( div [ E ) ) = — { g 0 e 0 


dt 


dE 

Ht 


Comme en dehors des charges 


on obtient finalement : 


div [ E ) = 0 


^ d 2 Ê > 

AÊ — g 0 e 0 -^- = 0 


( 12 . 5 ) 

( 12 . 6 ) 

( 12 . 7 ) 


Pour établir l’équat-ion aux dérivées partielles pour le champ magnétique B , calculons le 
rotationnel de chacun des membres du théorème d’Ampère-Maxwell : 


rot ^ rot = grad (div — AB = rot < g 0 e 0 


dE 

Ht 


mais 

div = 0 

et en inversant l’ordre des dérivations : 

— > d r— > 

-Aü = +^0^ rot [E 


or 


Donc 


ou encore 


rot ( E ) = 


dB 


-AB = -fi 0 e o 


dt 

d 2 B 

Ht 2 " 


s, q 2 ^ 

AB - h q eq = 0 

On obtient la même expression que pour le champ électrique. 
L’équation 


E 

B 


] = 0 

c 2 dE 1 B 


( 12 . 8 ) 

( 12 . 9 ) 

( 12 . 10 ) 

( 12 . 11 ) 

( 12 . 12 ) 

( 12 . 13 ) 


( 12 . 14 ) 


c = 


constitue l’équat-ion de propagation du champ électromagnétique dans le vide, où on a posé 
1 

V% £ o 
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12.2 L’onde plane progressive sinusoïdale 

12.2.1 Relation de dispersion 

L’onde plane progressive sinusoïdale est définie, en notation complexe, par : 

È(r,t) = Ê (12.15) 

où k est le vecteur d’onde donnant la direction de propagation de l’onde plane. 

En utilisant la défintion du laplacien vectoriel dans un système de coordonnées cartésiennes, 
on peut montrer que : 

A Ê = -k 2 Ê (12.16) 

L’équation de propagation s’écrit alors sous laforme 

a; 2 ! 

-/c 2 + — Ê = 0 (12.17) 

c z 

L’onde plane progressive sinusoïdale constitue une solution particulière de l’équat-ion d’onde 
seulement si la relation suivante, dite relation de dispersion, est satisfaite : 

k = ^ (12.18) 

12.2.2 Structure de l’onde uniforme plane 

L’onde plane progressive sinusoïdale doit également satisfaire le théorème de Gauss. On 
montre aisément que pour une onde plane progressive sinusoïdale : 

cliv = 0 est équivalente à — i k ■ Ê = 0 (12.19) 

Soit encore k-E = 0 ; ce qui revient à dire que le champ électrique E est perpendiculaire à la 
direction de propagation donnée par le vecteur d’onde k. Le champ électrique est dit transversal. 



y 
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12.3 Onde plane uniforme progressive et sinusoïdale 


12.3.1 Onde de polarisation rectiligne 

Jusqu’ici nous avons considéré un champ électrique E dont la variation en fonction de z 
et de t était quelconque. Nous allons maintenant examiner le cas où le champ E garde une 
direction constante (polarisation rectiligne) et varie sinusoïdalement ; il s’écrit alors : 


Ê 


Eq cos 


U 



E 0 cos luit — k • r 


où 

E 0 = vecteur constant 
u = pulsation de la fonction sinusoïdale 
k = lu/ c : constante appelée nombre d’onde 
k = k n : vecteur d’onde perpendiculaire au plan d’onde 

eut — k ■ f : phase instantanée ou plus simplement phase de la grandeur variable, 
c = - j== vitesse de propagation dans le vide . 

Les vecteurs forment toujours un trièdre direct et 

-> ~n x E k x E 

B = = 

c UJ 

Pour préciser cette onde, supposons qu’elle se propage suivant z' z d’où 

E = E 0 cos (c ut — kz ) 


On constate une double périodicité : 


Une périodicité temporelle : pour z donné le champ varie sinusoïdalement en fonction du 
temps avec une période 

T — — 


C jJ 


ou une fréquence 



eu 

27T 


(/ est en hertz). 

-Une périodicité spatiale : à un instant t donné le champ varie sinusoïdalement en fonction 
de z avec une période 

2tt _ A 
~~k~k 


(À est appelée la longueur d’onde dans le vide). On peut remarquer que la longueur d’onde 
À est égale à la distance parcourue par l’onde pendant une période. 
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12.3.2 Onde de polarisation quelconque 

Dans le paragraphe précédent, nous avons supposé que le champ E (donc B également) gar- 
dait une direction constante. Dans le cas général, il n’en est pas toujours ainsi et les composantes 
du champ peuvent se mettre sous la forme : 


E x = E 0x cos (eut — kz — 44) 

Ey = E 0 y COS (tut — kz — 44) 

0 i, 0 2 pouvant être différentes. 

Etudions le comportement du champ E dans le plan z = 0. Les résultats obtenus se retrouvent 
avec un décalage temporel dans tout plan z = cte. Les composantes du champ s’écrivent : 


E x = E 0x cos (tut - 44) 

Ey = E 0 y COS (tut — 44) 

E z = 0 

et si l’on prend pour origine des temps un instant où E x passe par sa valeur maximale on 
a : 


E x 

Eq x 


E 0 y 


cos (tut) 
cos (tut — 4) 


avec 4 — 02 0i - 

On peut déjà dire que l’extrémité du vecteur décrit une courbe inscrite dans un rectangle 

E 

de côtés 2E 0x et 2E 0y . D’autre part en développant l’expression de et en éliminant le temps 
il vient : 


Eo. 


E t 


0 y 


= cos (tut) cos (0) + sin (tut) sin (0) 


E ,, 


Ei 


0 y 


E X COS (0) + \ll - { ' E 


Ei 


Ox 


E, 


sin i 


y , 


r e e 

y 7 cos (0) 

2 


_E 0y E 0x 


\Eq x J 


sur 


E x 


Ei 


Ox 


En. 


Ei 


0 y 


E x E v 

2 f^f^ cos{ 

^0x y 


= sin 2 (0) 


Pour 0 quelconque, cette équation est celle d’une ellipse : on dit que l’onde a une polarisation 
elliptique ; pour 0 = mn(m = 0, 1, 2, ... ) l’ellipse dégénère en une droite et l’onde est dite à 
rectiligne. Enfin si E 0x = E 0y et si 0 = (2/ + l)n/2 l’onde est dite à polarisation circulaire. 
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12.3.3 Notation complexe 

Considérons une onde plane sinusoïdale polarisée rectilignement se propageant suivant une 
direction n et ayant pour vecteur d’onde k = k n . 



Le support du vecteur champ électrique E de direction constante a pour vecteur unitaire u 
et nous pouvons écrire : 

Ê = Ë 0 cos (ut — k ■ r + ip'j 

sous la forme 

Ë = E 0 cos (ut — k • r + ip] u 


Avec Eq = cte =amplitude du champ E. 

On sait qu’une grandeur sinusoïdale peut être représentée par un nombre complexe et qu’à 
chaque instant cette grandeur est obtenue en prenant la partie réelle du complexe qui la repré- 
sente. Ainsi dans le cas présent, nous pouvons représenter la mesure 

E = En cos (eut — k ■ r + 


du vecteur champ électrique sur son support par le nombre complexe 

E* = Eq exp (ut — k ■ r + 'ip) 
et nous avons bien à chaque instant 

E = Re (E*) 

Ce complexe appelé la composante complexe de E sur son support peut aussi bien s’écrire : 

E* = En e Ë- tp+ ^)e jut = Ë e jujt 


E = En e j (- lF +Ë 
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est l’amplitude complexe de la composante de E sur a . Finalement nous caractérisons une 
onde plane polarisée rectilignement par sa représentation complexe : 


Ê* = E 0 e j ^ 
B* = B 0 e ji> 


12.4 Energie électromagnétique : vecteur de Poynting 

La propagation de l’énergie se manifeste expérimentalement dans de nombreux cas : 

— On peut ressentir son effet si l’on s’expose aux rayons solaires ou au rayonnement d’une 
source chaude ; 

- De même tout émetteur radio expédie de l’énergie à travers l’espace, une infime partie de 
cette dernière étant captée par votre récepteur radio. 

Nous allons essayer de relier localement cette énergie qui se propage, au champ électro- 
magnétique qui la transporte. Nous supposerons le milieu de propagation parfait, c’est à dire 
homogène, isotrope et linéaire. 


12.4.1 Onde de forme spatiale et temporelle quelconques 

Nous admettrons que les densités d’énergie électrique et magnétique calculées en régime 
stationnaire sont toujours valables en régime variable ; la densité d’énergie électromagnétique 
w en un point quelconque du milieu parcouru par une onde électromagnétique est donc à chaque 
instant : 

w = densité d’énergie électrique + densité d’énergie magnétique 




Considérons dans le milieu, un volume r limité par une surface ( S ). L’énergie électroma- 
gnétique qu’il contient est à chaque instant : 


W = 



Pendant un temps dt l’accroissement d’énergie dans (r) sera dW et la puissance instantanée 
p' acquise par ce volume sera 

, dW 


P = 


dt 



dw , 

m dT 


( T ) 


On a : 


\ dB 

^ 

B 

' “ dt 

et rot 

ho 


— £o' 


dE 

dt 
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donc 


dw 

Ht 


= E ■ 


rot 


B 

Mo 


B 

IM) 


■ rot E 


D’après une relation de transformation, on a : 

div 


donc 


et 


Ex — 
Mo 


= — • rot (e) - E ■ rot ( — 


Mo 


Mo 


diu / - B 

—— = —div \ E x — 

ât \ /i 0 


P = 



(t) 


B 

div | Ë x — | c/r 
Mo 


La puissance électromagnétique instantanée perdue par le volume (r) est : 


~P = 



(t) 


B 

div | Ë x — | c/r 
Mo 


Elle représente la puissance électromagnétique qui sort du volume (r), c’est à dire la puis- 
sance moyenne p rayonnée par ce volume. 


P = 



(t) 


B 

div | Ex — | c/r 
Mo 


D’après la formule d’Ostrogradsky, on peut écrire : 


Le vecteur 



• dS = 



P-dS 


(S) 


R 


Ê 


B 
x — 
Mo 


est appelé le vecteur de Poynting. Sa direction donne en chaque point, la direction d’écoule- 
ment de l’énergie et son flux à travers une surface est égal à la puissance électromagnétique 
instantanée rayonnée par cette surface. Les courbes tangentes en chaque point au vecteur de 
Poynting peuvent être considérées comme des trajectoires de l’énergie ; on les appelle les ‘hayons 
électromagnétiques” . 
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12.4.2 Onde plane progressive et uniforme sinusoïdale 

Puisque (È, B, nj forment un trièdre trirectangle direct le vecteur R a même direction et 

sens que k c’est à dire que l’énergie s’écoule dans le sens de propagation (ce résultat n’est pas 
général ; en effet dans un milieu anisotrope par exemple R et k ne sont pas colinéaires) . 

La puissance instantanée p„ traversant une surface unitaire ( S ) perpendiculaire à la direction 
de propagation est 


Pu= J J R ■ ds = 

( 5 ) ( 5 ) 


R 


dS = 


R 


clS = 


R 


S 


( 5 ) 


La puissance moyenne traversant (S) est alors 


or B _L É, et 


B 


T 


— \Jdo £ 0 


i r T 


1 

f T 

- B 

dt = — 

R 

dt = — 


Ex — 

T Jo 


T 

Jo 

ho 


dt 


E 


, d’où 


(Pu} = 


1 

T 


/ E\^dt 

I o V / J o 


Si l’onde est polarisée rectilignement alors 

E = E 0 cos [ut — k ■ r 


{Pu) 


TP f \ o cos2 (ut — k-r\ dt 

T Jo V /‘o V J 


(Pu) 



f T 1 

r 

2 (ut — k ■ J ) 

/ h 

1 + COS 

Jo 2 

- 

. V / J . 


dt 




PE 


où E e ff= valeur efficace de E. 

Le flux d’énergie traversant par unité de temps l’unité de surface perpendiculaire à la di- 
rection de propagation est une constante dépendant du milieu et proportionnelle au carré de la 
valeur efficace du champ électrique. 
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Chapitre 13 


Réflexion et transmission des ondes 
élect r omagnét iques 


Nous allons dans ce chapitre utiliser les équations de maxwell pour étudier la propaga- 
tion d’une onde d’une onde électromagnétique dans un milieu isolant puis sa réflexion et sa 
tranmission à un interface. 


13.1 Equations de Maxwell dans les milieux parfaits 

Nous supposerons que l’onde se propage dans un milieu illimité satisfaisant les conditions 
suivantes ( milieu parfait) : 

•Le milieu est considéré comme continu ; ce qui est légitime pour les ondes habituellement ren- 
contrées (hertziennes et lumineuses) dont les longueurs d’onde sont, en général, très supérieures 
aux distances inter-atomiques. 

•Le milieu est homogène et isotrope 

•Le milieu peut être caractérisé par une permittivité £ et un permabilité /i et une conductivité 
7 constantes égales à celles définies en régime stationnaire (milieu linéaire). 

Ces hypothèse simplificatrices permettent de traiter un grand nombre de problèmes ; il faut 
toutefois signaler qu’elles ne tiennent pas compte des phénomènes d’absorption, d’hystérésis et 
de dispersion ( £ , // et 7 peuvent alors être complexes et dépendre de la fréquence). 

Dans le cas des milieux homogènes, linéaires et isotopes, les équations de Maxwell de- 
viennent : 

- Théorème de Gauss pour E 


Forme locale 

Forme intégrale 

div (Ê) = Plibres 

J 

f J Ê ■ dS = j If f>libres dr 

(S fermée) (r) 


- Théorème de Gauss pour B 
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Forme locale 

Forme intégrale 

div (^b) 

1 =0 

J 

f J B ■ dS = 0 

( S fermée) 


- Loi de Faraday 


Forme locale 

Forme intégrale 

rôt 

i)B 
1 “ ~~dt 

J 

(C 

( Ê-cü = -j t JJ^-dS 

) (S) 


- Théorème d’Ampère-Maxwell 


Forme locale 

Forme intégrale 

rot (^Ê'j 

dÈ 

1 = h 3 + PZ-gjT 

j 

(r 

lë-dî= JJ | 

) (5) 

( ■ BÉ) 

^ 3 + f,E SF) 

I • dS 


où 

- £ = e 0 £ r est la permittivité absolue du milieu diélectrique 

- £ r est la permittivité relative du milieu diélectrique. 

- £q = 8.854187817 x 10~ 12 F m” 1 est la permittivité du vide 
/i = Pohr es t la perméabilité absolue du milieu 

- /i r est la perméabilité relative du milieu. 

- /i 0 = 1.2566370614 x 10~ 6 NA~ 2 est la perméabilité magnétique du vide 


En tout point d’un conducteur, il existe une relation entre le vecteur densité de courant j et 
le champ électrique total E ( somme du champ électrostatique Es et du champ électromoteur 
E rn ) dite relation d’Ohm Kirchhoff : 


j = 7 (Es + 

3 = 7 È 

où 7 est la conductivité du milieu conducteur. 

Cas particuliers : 

- Isolants : 7 = 0 

- Conducteurs parfaits : 7 — > 00 

13.2 Propagation dans les milieux diélectriques 

Les milieux diélectriques sont des mieux isolants. Leur conductivité est extrêmement faible, 
de l’ordre de 10 20 à 10 12 (12 -m -1 ), celle d’un conducteur métallique étant de l’ordre de 
10 7 (12 • mT 1 ) , à température ambiante. Il est donc tout à fait raisonnable de prendre pour 7 
la valeur 7 = 0. Par ailleurs dans un tel milieu, Pi ibre = 0. 

Les équations de Maxwell se simplifent alors en : 
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div 

\ É ) 

= 0 

div 

(5) 

= 0 



(ê) 

dB 

rot 


V J 

dt 

> 

(5) 

dÊ 

rot 

= » £ ~dt 


En utilisant la même démarche que dans le chapitre précédent, on peut montrer que le 
champ électrique et le champ magnétique satisfont les équations de propagation suivantes : 


A E 


A B 


1 d 2 Ê 
V 2 'dt 2 
1 d 2 B 


V 2 dt 2 

où la vitesse de propagation de l’onde est : 


0 

0 


yj . \J V 0 £ 0y/ f^r^r Tl 

n = v /// r c r est l’indice de réfraction (ou indice optique) du milieu. Dans les milieux réels n 
est constant pour les grandes longueurs d’onde, tandis que pour les hautes fréquences il faut 
faire intervenir le phénomène de dispersion qui entraîne une dépendance de n avec la fréquence. 
Dans la plupart des diélectriques n r = 1, d’où n = sfE~ r . 

On peut égalment montrer que l’impédance caractéristique d’un tel milieu peut s’écrire : 


Z = 

où Z 0 est l’mpédance caractéristique du vide. 


Zo 

n 


13.3 Relations de passage 

Composante tangentielle et composante normale de E 

A la traversée d’une surface (E) séparant deux milieux et portant des charges vraies avec 
une densité supercielle a, les relations locales s’écrivent : 


Et2 — Eti 
£2 En 2 £1 Ejsri 


0 

a 


où E t est la composante de E dans le plan tangent à (E) en M , tandis que E N mesure de 
la composante de E suivant la normale n en M et orientée du milieu (1) vers le milieu (2). 
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Composante tangentielle et composante normale de B 

A la traversée d’une surface (S) séparant deux milieux et parcourue par des courants vrais 
de densité superficielle js , les relations locales B deviennent : 

Bn2 = B N i 

Ih Bt 2 ~ ^ i Bti = js x n 

où B n est la composante de B suivant la normale au point M considéré orientée du milieu 
(1) vers le milieu (2);tandis que B T est la composante de B dans le plan tangent à (S) en M. 

13.4 Lois de Snell-Descartes 

Considérons deux milieux 1 et 2, parfaits, isolants, non magnétiques et semi-infinis, carac- 
térisés respectivement par S\ ,/.i 1 = /j 0 , = 0 et par £ 2 ,h 2 = l L o> 72 = et séparés par une 

interface plane. On constate expérimentalement qu’une onde plane incidente arrivant à cette 
interface donne généralement naissance à deux ondes, une onde plane réfléchie et une onde 
plane transmise ou réfractée ( la réflexion entre deux milieux parfaits et isolants est souvent 
appelée réflexion vitreuse ) 

Nous nous proposons de déterminer entièrement les ondes réfléchie et transmise en appli- 
quant les équations de propagation et en tenant compte des conditions aux limites. Quel que 
soit l’état de la polarisation de l’onde incidente, il est toujours possible de la décomposer en 
deux ondes à polarisation rectiligne perpendiculaires entre elles en projetant le vecteur inci- 
dent sur deux axes perpendiculaires ; nous étudierons donc d’abord le cas d’une onde incidente 
polarisée rectilignement. 

Définissons un trièdre de référence par une origine O appartenant au plan séparant les deux 
milieux, l’axe Oz étant perpendiculaire à ce plan et l’axe Oy étant tel que le plan yOz contienne 
le vecteur k\ de l’onde incidente. En appelant u 1 , u \ et u 2 les vecteurs unitaires des supports 
des vecteurs champs électriques des différentes ondes, nous pouvons représenter ces dernières 
en notation complexe par 

- Onde incidente : Ei = E 0Ï Tt 1 k v r 0 

- Onde réfléchie : E[ = E’ m 

- Onde transmise : E 2 = E 02 ~u 2 e J ( a;2t ~ A:2 ' r2 ) 

Dans ces expressions E m , E' 01 et E 02 sont les amplitudes complexes des différents champs 
électriques. En un point quelconque M 0 du plan xOy tel que OM 0 = r o,les ondes sont 
représentées par 

- Onde incidente : Ei = E m ~u 1 e i( Wi *- k i- r 0 ) 

- Onde réfléchie : E[ = E’ m d?j 

- Onde transmise : E 2 = E 02 IÎ 2 e J ( a;2t_A:2 ' ro ) 
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Réflexion-transmission 


La continuité de la composante tangentielle du champ entraîne une relation entre les com- 
posantes tangentielles de la forme 

~a 1 E 01 A Ult - tl - r °) + EoxA^-^ 0 ) = ~a 2 

Pour que cette relation soit satisfaite à chaque instant en tout point M 0 de la surface de 
séparation, il faut que 


c d\t — k i • r o = L 0 \t — k [ ■ r o = co 2 1 — k 2 ■ r 


ce qui implique 


uji = u4 = cu 2 et que k i • r 0 = k \ ■ r 0 = k 2 ■ r 0 
La dernière relation conduit à 

( k[ - k i) • rqj = 0 

(k 2 — k -~r o = 0 

qui seront satisfaites quel que soit le point M 0 si ( ^ k \ — fcxj et (^k 2 — k l'j sont per- 
pendiculaires au plan xOy, c’est à dire si k ±, k \ et k 2 sont coplanaires : ces trois vecteurs 
appartiennent au plan d’incidence ( k i,’e’ 3 V 
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En prenant r 0 = e 2 , on obtient 

k i • ~e 2 =t- k[ sin (i\ ) = k\ sin (ti) 
k i ■ ~e 2 => k2 sin (i 2 ) = k\ sin (ti) 

UJ LU UJ 

or ki = — , k[ = — et k 2 = ; d où 

Vl Vl V 2 

sin (ij) = sin (ii) =>• i[ = ii 
sin (i 2 ) sin (ii) 

V 2 = ^ 

En multipliant les deux membres de la dernière relation par c (vitesse de propagation dans 
le vide) et en faisant intervenir les indices des milieux, on obtient 

ni sin (i 2 ) = n 2 sin (ii) 

En définitive nous avons montré que : 

- Les ondes réfléchie et transmise ont même fréquence que l’onde incidente 

- Les vecteurs d’onde réfléchi et transmis sont dans le plan d’incidence (plan et la normale 
à la surface de séparation au point d’incidence) 

- L’angle de réflexion i\ est égal à l’angle d’incidence i\ 

- Les angles de transmission i 2 et d’incidence fl sont liés par ni sin (t 2 ) = n 2 sin (ti) (loi de 
la réfraction). 

Ces résultats représentent les lois de Snell-Descartes pour la réflexion et la réfraction. 

13.5 Formules de Fresnel 

Nous allons étudier successivement le cas d’une onde dont le champ électrique est polarisé 
rectilignement dans le plan d’incidence puis le cas d’une onde dont le champ électrique est 
polarisé perpendiculairement au plan d’incidence 

13.5.1 Champ électrique dans le plan d’incidence 

Le trièdre de référence est choisi comme précédemment et les champs positifs des champs 
électrique et magnétique sont indiqués sur la figure (Pour satisfaire les raisons de continuité 
et de symétrie, il faut supposer que les ondes réfléchie et transmise ont leur vecteur champ 
électrique polarisé dans le plan d’incidence). 

Sachant que les différents vecteurs unitaires u \. u \ et u 2 sont dans le plan d’incidence, les 
différentes ondes ont pour représentation complexe : 

- Onde incidente : Efl = E m ~u 1 e K UJlt ~^ 1 " r x ) 

- Onde réfléchie : E[ = E' m dfl; e i("P-^i) 


k i • "e’ 2 = 
k 2 -~e 2 = 
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- Onde transmise : E 2 = E 02 ~U2 e j ( U2t k2 ' r 2 ) 



Champ électrique parallèle au plan d’incidence 


Dans le milieu l l’onde résultante est la somme de l’onde incidente et de l’onde réfléchie. Les 
continuités des composantes tangentielles du champ électrique et de l’excitation magnétique à 
la traversée de la surface séparant les deux milieux permet, si l’on tient compte de // , = /j 2 = // 0 . 
de B 0 = Eq/V et de la condition k y r 0 = k [ ■ r 0 = k 2 • r 0 d’obtenir les relations 

E 0 i cos (ii) + E ' 01 cos (A) = E 02 cos (i 2 ) 

-E'01 _ E ' 01 _ Eq2 

“ 1 ¥ 

(le signe moins provient du sens de B\ opposé à celui de Ox ; ce sens a été choisi ainsi pour 
respecter le fait que (^Ë. B , k \ est toujours direct). 

A partir des relations ci-dessus, on obtient E[ n et Eq 2 en fonction de E (n : 


or 


-E'01 


cos (A) , 
cos(i 1 ) + 01 


E, 


02 


Vy 

Vo 


E' 


01 


^02 


cos ( i 2 ) 
cos (A) 


-E'01 


Vi _ n 2 _ sin(A) 
V 2 ni sin ( i 2 ) 
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d’où les formules de Fresnel 

2 cos (ii) sin (i 2 ) 2 cos (ii) sin (i 2 ) 

sin (ii) cos (ii) + sin (i 2 ) cos (i 2 ) sin (ii + i 2 ) cos (ii — i 2 ) 

sin (ii) cos (ii) — sin (i 2 ) cos (i 2 ) tan (ii — i 2 ) 

sin (ii) cos (ii) + sin (i 2 ) cos (i 2 ) tan (i x + i 2 ) 

Nous notons [ ]y pour indiquer que le champ électrique est dans le plan d’incidence. 

13.5.2 Champ électrique perpendiculaire au plan d’incidence : 

Les dispositions et les sens positifs pour les champs sont indiqués sur la figure cidessous 




Champ électrique perpendiculaire au plan d’incidence 

si l’on tient compte de /j , = /i 2 = /i 0 . de B 0 = E 0 /V et de la condition /c i • = k' 1 -~r 0 = 

k 2 ■ r o les conditions de continuité à la traversée de la surface séparant les deux milieux 
permettent d’obtenir les relations 

E 01 + E' m = E 02 

— — cos(îi) + — cos pi) = — — cos(z 2 ) 
vi V\ V 2 

En tenant compte 

Vj_ _ ^2 
V 2 n i 

on obtient finalement les formules de Fresnel 


sin (L) 
sin (i 2 ) 
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T ± 


R± 


Eq2 


Eqi 



Eqi 


2 sin (i 2 ) cos (ii) 
sin (ii + i 2 ) 
sin (ii - i 2 ) 
sin (ii + i 2 ) 


Nous notons [] ± pour indiquer que le champ est perpendiculaire au plan d’incidence 


13.5.3 Discussion des résultats 

Dans les formules de Fresnel, les coefficients de réflexion et de transmission sont des nombres 
réels ; il en résulte que le déphasage introduit par la réflexion ou la transmission aura toujours 
une valeur nulle ou égale à n. 


Incidence normale : 


Nous obtiendrons les résulats de l’incidence normale en faisant tendre i; vers zéro et en 
tenant compte alors de ni i\ = n 2 i 2 . Les relations de Fresnel s’écrivent alors 


T\ 

Ri 


T\ 


R± 


2 ni 

ni + n 2 
ni - n 2 
ni + n 2 


- Les ondes réfléchie et transmise sont indépendantes de la polarisation de l’onde incidente 

- L’onde transmise est toujours en phase avec l’onde incidente 

- La composante électrique est déphasée de n à la réflexion si ni < n 2 ; elle n’est pas 
déphasée si ni > n 2 . 

Intéressons nous maintenant au comportement énergétique de l’onde à l’interface plane ; ce 
comportement sera caractérisé par les deux coefficients suivants 

T . puissance réfléchie 

- Le pouvoir réflecteur an = 

puissance incidente 

. . puissance transmise 

- Le facteur de transmission «t = 

puissance incidente 

Les relations donnant les intensités, c’est à dire les puissances par unité de surface, trans- 
portées par les différentes ondes sont données par 


_ £pi»-l . rf 

2 en ’ 1 


E&ni 

2c/h 


et I 2 


Ep 2 n 2 
2 c/x 2 


avec Un 


t-h — ho 


d’où 


OiR 


a x 


'Ki 

_Eqi 

Eq 2 

Eqi 


“1 2 


“1 2 


ni -n 2 1 2 
ni + n 2 _ 

4 nin 2 

(ni + n 2 ) 2 
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Remarque : La conservation de la puissance implique Or + cpr = 1 ce que l’on peut vérifier 
à partir des deux dernières relations. D’autre part, le pouvoir réflecteur et le facteur de trans- 
mission ne dépendent pas du sens de propagation de l’onde incidente (milieu 1 — > milieu 2 ou 
milieu 2 — > milieu 1). Pour les ondes lumineuse, dans le cas où les milieux sont l’air et le verre 
(ni = /, ri 2 = 1.5) on a or = 0.04 et «t = 0.96. 

Incidence oblique 

Nous allons discuter le comportement de l’onde en fonction de son angle d’incidence q 
(0 < ii < |) dans le cas particulier où ni < n 2 : 

- Onde transmise : Puisque sin(i 2 ) = ^sin(p) on a î 2 < i\ et il existera toujours une onde 
transmise. Il en résulte que le champ électrique de l’onde transmise, aussi bien dans le 
cas de polarisation dans le plan d’incidence (noté ||) que de polarisation perpendiculaire 
au plan d’incidence (noté _L), reste en phase avec celui de l’onde incidente sur la surface 
de séparation 

- Onde réfléchie : Pour une onde polarisée perpendiculairement au plan d’incidence (_!_) le 
champ électrique réfléchi est toujours en opposition de phase avec celui de l’onde incidente. 
Pour une onde polarisée dans le plan d’incidence (||) et sous incidence normale, le champ 
électrique réfléchi est en opposition de phase avec le champ incident. Lorsque l’angle 
d’incidence atteint la valeur i\B définie par %\b + *2 = § , le champ électrique réfléchi 
s’annule. Cet angle %\b est appelé angle de Brewster et en tenant compte de «/ sin(q) = 
n 2 sin(i 2 ), on trouve qu’il a pour valeur %\b = atan(n 2 /?p). Pour un angle d’incidence 
proche de | (incidence rasante) ; le champ électrique réfléchi a même amplitude que le 
champ incident et ils sont en opposition de phase. 

Si l’onde lumineuse arrive directement d’une source (le soleil par exemple) sur une surface 
vitreuse, elle est polarisée elliptiquement et on peut la considérer comme la superposition 
de deux ondes polarisées rectilignement l’une dans le plan d’incidence l’autre perpendicu- 
lairement au plan d’incidence. Suivant l’angle d’incidence, la lumière réfléchie sera plus ou 
moins polarisées rectilignement cette polarisation étant totale pour l’angle de Brewster 
la vibration lumineuse est alors rectiligne et parallèle à la surface réfléchissante. 

- Pouvoirs réflecteurs : Les pouvoirs réflecteurs (intensité réfléchie sur intensité incidente) 
ont respectivement pour valeur 


«iî|| 

®R± 


tan 2 (ii — i 2 ) 
tan 2 (ii + i 2 ) 
sin 2 (ii - i 2 ) 
sin 2 (ii + i 2 ) 


La figure dessous donne pour une onde lumineuse et dans le cas air- verre, la variation des 
pouvoirs réflecteurs en fonction de l’angle d’incidence (i\ b = 56,3°) . 


http: / / dj elouah.ifrance.com 





13.6 Réflexion sur un conducteur parfait 


109 



Angle d'incidence 


Variation du pouvoir réflecteur en fonction de l’angle d’incidence 


13.6 Réflexion sur un conducteur parfait 


Considérons une onde plane incidente uniforme sinusoïdale polarisée rectilignement ( E / || 
Ox) se propageant dans le vide ( ou l’air) suivant zO et arrivant sous incidence normale, à la 
surface plane d’un conducteur de conductivité infinie (conducteur parfait). Des considérations 

énergétiques montrent qu’il n’y pas d’onde transmise i^E 2 = 0 etB 2 = 0 j et que les seuls 
courants vrais pouvant être induits par cette onde le sont sur la surface du métal. Les raisons 
de symétrie impliquent que les directions de polarisation des vecteurs incident et réfléchi sont 
identiques. 



Réflexion métallique en incidence normale 
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Les champs électrique et magnétique des ondes incidente et réfléchie sont 

~E ! = E 01 ~e 1 ]l3 l = e j{ “ t+kiz) ~e 2 

c 

E[ = E' m e j{wt - kiz) ~e 1 -,^ 1 = ^ e j{ut - kiz) ~e 2 

car ~u i = ~u\ = ~6i , ~v' \ = d?', =e2; /c ! = — &u ~ë *3 ; = ~r\ = 2 ^3 et fî 0 = —• 

c 

La continuité de la composante tangentielle du champ électrique permet d’écrire à la tra- 
versée de la surface de séparation (z = O) 

{E 01 + E' (n } z=Q = 0 d’où Eqi = —E' 0 1 


Le champ électrique réfléchi a même amplitude que le champ électrique incident et il est 
déphasé de n par rapport à ce dernier. Pour satisfaire à la directivité du trièdre (^E , B , k j le 

champ magnétique est réfléchi suivant B 1 c’est à dire que sa réflexion s’effectue sans changement 
de phase. En définitive les ondes incidente et réfléchie s’écrivent 


Et = E m e j{ “ t+kiz) ~et ; ~Bt = e j(ujt+kiz) ~e 2 

c 

E\ = —E 0 i e j{u}t ~ kiz) x ; ~Bt = e ~£ 2 

c 

Dans le vide, l’onde résultante est la somme de l’onde incidente et de l’onde réfléchie et ses 
vecteurs champs ont pour valeur 

E = Et + ~E[ = E 01 ( e +jklZ - e ~ jklZ ) e jut ~et 
E = 2E 01 sm(ktz) e^+ï) ~t 1 


~B = Bt + B't = (e +jkiz + e~ jkiz ) e jljjt 

c y ' 

~B = 2^ cos (ktz)e j ^ t+n) ~et 

En revenant à la notation réelle 

E = 2 E 01 sin (ktz) cos (uit + — j ~e 1 
^ 

B = 2^— cos (ktz) cos (ut + tt) et 

Ces relations montrent que les champs ne se propagent plus mais qu’ils oscillent sinusoïdale- 
ment en fonction du temps avec une amplitude qui est fonction de la distance z, leur déphasage 

7 T 

étant de — on dit que 1’ onde est stationnaire. 

Z 
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L’amplitude du champ électrique est nulle pour k\Z = ‘=fz = l i r(7 = 0,1,2, ) , c’est à 

dire pour z = 

L’amplitude du champ magnétique est nulle pour k\z = ^ z = (l + |) n (1 = 0, 1, 2, ) , 

c’est à dire pour 2 = (/ + |) |. 

Les points où l’amplit-ude est nulle sont appelés les noeuds et les points où l’amplitude est 
maximale sont appelés les ventres. 
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Equations différentielles 


A.l Introduction 

Les oscillations linéaires des systèmes à un degré de liberté sont régies par des équations 
différentielles du second ordre à coefficients constants. Une équation différentielle linéaire du 
second ordre à coefficients constants est une relation entre une variable dépendante y et ses 
dérivées première et seconde par rapport à une variable indépendante t, qui peut s’écrire sous 
la forme : 


+ 2 5^-+ou 2 y = A(t) 


d 2 V , nx dy 
dt 2 + dt 

Les coefficients <5 et cuo sont des constantes réelles positives. Dans les problèmes de vibration 
le paramètre t représente le temps et on note par convention : 


cly 

dt 

fy 

dt 2 

D’où l’écriture de l’équation différentielle : 


ÿ 

ÿ 


ÿ + 25ÿ + uj 2 0 y = A(t) 


A. 2 Equation homogène 

L’équation différentielle est dite homogène si A(t) = 0 : 

ÿ + 2 5 ÿ + u 2 0 y = 0 

On recherche des solutions qui sont de la forme y(t) = e st . Dans ce cas l’équation différen- 
tielle s’écrit : 
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[s 2 + 28 s + cüq] e st = 0 

Cette équation doit être vraie quel que soit t. ce qui implique : 

s 2 4- 2 S s cUg = 0 

On obtient ainsi une équation du second degré en s dite équation caractéristique. Les racines 
de cette équation caractéristique sont : 

Si = — 5 + \J cuq — S 2 

s 2 = -6 - \JujI~8 2 

La solution de l’équation différentielle s’écrit alors : 

y(t ) = Aie sit + A 2 e S2t 

Ai et A 2 sont deux constantes d’intégration que l’on peut déterminer à partir des conditions 
initiales : 


yit = 0 ) = 2/0 

ÿ{t = 0 ) = ÿo 

Sachant que 8 et co 0 sont des nombres réels positifs, Si et s 2 sont négatives ou complexes 
avec une partie réelle négative. La nature des solutions si et s 2 de l’équation caractéristique 
dépend de la valeur relative de S par rapport à cuo . Ainsi trois cas sont à envisager : 

A. 2.1 Régime fortement amorti ( ô > Uq ) 

Dans ce cas, les deux racines si et s 2 sont réelles et négatives. L’écriture des conditions 
initiales donne un système de deux équations : 


Ai + A 2 — yo 
siAi + s 2 A 2 = ÿo 

dont les solutions sont les constantes d’intégration Ai et A 2 : 


Ai 

A2 


S2Vo - ÿo 

S2 - Si 

ÿo - Siyo 

S2 ~ Si 
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y(t) est une fonction décroissant sans oscillations vers zéro lorsque t augmente. La forme 
exacte de y(t) dépend des valeurs de Ai et A 2 qui sont déterminées par les conditions initiales. 
Pour les conditions initiales particulières suivantes :yo ^ 0; i/o = 0, la figure suivante représente 
les variations de y au cours du temps 



A. 2. 2 Régime critique ( <5 = c uo ) 

L’équation caractéristique possède une racine double : 

si = s 2 = -S 

La solution générale de l’équation différentielle est de la forme : 

y(t ) = (Ai + A 2 t) e St 

Les constantes d’intégration Ai et A 2 sont déterminées à partir des conditions initiales et 
valent : 


Ai = y 0 

A2 = ÿo + à Vo 

y(t ) et une fonction décroissant vers zéro quand t augmente. Il est aisé de vérifier que cette 
situation correspond à la décroissance la plus rapide de y(t). Ce cas correspond au régime 
critique. Pour le cas particulier de conditions initiales :y 0 7 ^ 0; y 0 = 0, la solution est : 

y{t) = Vo e~ st (1 + St) 

La figure ci-dessous représente les variations de y(t). 
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A. 2. 3 Régime pseudo-périodique ( ô < ujq ) 

Dans ce cas si et s 2 sont complexes conjugués : 

si = -5 + j^ul - S 2 
s 2 = -5-j^Juj l - 6 2 

où j représente le nombre imaginaire pur vérifiant la relation j 2 = —1. Posons : 

uj d = \JujI~ S 2 

La solution générale s’écrit alors : 

y(t) = e~ 5t [Ai e j “ Dt + A 2 e~ juJDt ] 

Sachant que : 

e ju D t _ cos (uj D t) + j sin (oj D t) 
e ~ju D t _ cos(uj ot) — j sm(uj ot) 

y(t ) s’écrit : 

y{t) = e~ 5t [(Ai + A 2 ) cos(uj D t) + j (Ai - A 2 ) sin (u D t)] 
y if) étant réelle, les nombres complexes Ai et A 2 doivent vérifier les relations : 
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A 1 + A 2 : réel 

Ai — A 2 : imaginaire 


A\ et A 2 sont donc complexes conjugués et peuvent se mettre sous la forme : 


Ai = A'e jlfi 
A 2 = A! e~ jip 


y(t ) s’écrit alors : 


y(t) = A'e~ st [ e A“ Dt +A + e -j^Dt+ v y 


soit finalement, en posant A = 2 A' : 


y(t) = A e St cos(c ont + <p) 

y(t ) peut être interprétée comme une fonction périodique de pulsation üjd . de phase initiale 
ip et d’amplit-ude Ae H décroissant exponentiellement au cours du temps. On peut définir une 
pseudo-période 


CÜD 

A et (p sont deux constantes d’intégration définies à partir des conditions initiales : 

y 0 = A cos((f) 

ÿo = —A cos (ip) [ô + oj d tan(</?)] 


D’où l’on tire : 


<P 

A 


— arctan 


up y o 
ÿo + 5yo 


x/(u D y 0 ) 2 + ( ÿ 0 + ôy 0 ) 2 


UJp 


La figure ci-dessous représente les variations de y(t) dans le cas particulier de conditions 
initiales :yo ^ 0; ÿo = 0. 
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Remarque : Cas particulier où 5 = 0 : 
L’équation différentielle s’écrit : 


ÿ + uj 2 0 y = 0 

La solution s’exprime dans ce cas : 

y(t) = A cos(to> 0 t + 4>) 

Cette solution est appelée solution harmonique car y(t) est une fonction sinusoïdale du 
temps, de pulsation üJ(, , de période T 0 = 2tt/oj 0 et dont l’amplitude A et la phase initiale à 
sont déterminées par les conditions initiales : 


y 0 = A cos(0) 
ÿo = -uJq A sin (0) 


d’où l’on tire : 


0 

A 


— arctan 


k'o Vo 

ÿo 




2 


La figure suivante représente les variations au cours du temps de y(t) dans le cas particulier 

où : i/o = 1 et ÿ 0 = 0. 
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A. 3 Equation avec second membre 

A. 3.1 Solution générale 

Soit y(t) la solution générale de l’équation différentielle avec second membre : 

ÿ + 2Sÿ + ujly = A(t ) 

Soit z/nit) la solution de l’équation homogène et soit 'yp(t) une solution particulière de 
l’équation avec second membre ; y h A) et yp(t) sont les solutions respectives des deux équations 
différentielles suivantes : 


ÿ)H + 2 S ijH + wIvh = 0 
ÿp + 2 S ÿp + cuq yp = A(t ) 

Les opérations de dérivation qui interviennent étant des opérations linéaires, l’addition 
membre membre des deux équations différentielles précédentes donne : 


(Vh + ÿp) + 2 S ( ÿ H + ÿ P ) + ujI (y h + y P ) = A(t ) 

Ainsi la solution générale peut être obtenue en faisant la somme de la solution homogène 
et d’une solution particulière : 


vit) = y H (t) + yp{t) 
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A. 3. 2 Cas particulier où A(t) est constante 

Considérons le cas particulier important où A(t) = A 0 , A 0 étant une constante réelle. 
L’équation différentielle s’écrit alors : 


On peut vérifier que 


ÿ + 2 S ij + CÜQ y = A) 


^-0 
y = 77 

UJ 0 

constitue une solution particulière de l’équation avec second membre. Selon le cas la solution 
générale de l’équation avec second membre sera : 

- Régime fortement amorti (S > ouo) 


y(t) — Aie Slt + A 2 e S2t H — ^ 


ws 


Régime critique (S = cu 0 ) 


y(t) — (Ai + A 2 t) e 5t H j 

uj 0 

- Régime pseudo-périodique (S < coq) 

y(t ) = A e~ St cos(a j D t + ip) + ^ 

uj 0 

Dans les trois cas considérés les constantes d’intégration A 1 , A 2 , A et sont déterminées à 
partir des conditions initiales. Considérons à titre d’exemple le cas particulier y 0 = 0 et ÿo = 0 . 
La résolution des systèmes d’équations obtenus permet de calculer les constantes d’intégration 
et de tracer les variations de y(t) pour les trois régimes. Les variations correspondantes de y(t) 
sont représentées par la figure ci-dessous. 
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A. 3. 3 Cas particulier où A(t) = Aq cos (fit) : 

Nous pouvons vérifier que yp(t) = Y 0 cos (Qt + 0) constitue une solution particulière de 
l’équation différentielle avec second membre à condition que l’amplitude Y 0 et la phase 9 véri- 
fient la relation : 


La solution complète s’écrit alors suivant le cas : 

- Régime fortement amorti ( S > cuo) 

y(t) = Aie sit + A 2 e S2t + Y 0 cos (fili + 9) 

- Régime critique (5 = cno) 

y(t) = (Ai + A 2 t) e~ 5t + Y 0 cos (fit + 9) 

- Régime pseudo-périodique (<5 < a>o) 

y(t) = A e~ ôt cos(a j D t + ip) + Y 0 cos (fit + 9) 

Comme dans le cas précédent les constantes d’intégration A t . A 2 . A et 0 sont déterminées 
à partir des conditions initiales. La figure suivante représente le résultat obtenu dans le cas 
particulier où :( yo = 0 , ÿo = 0 ). 


[cüq — fil 2 ] Y 0 cos (Qt + 9) — 2 SQ sin (Qt + 9) = A 0 cos(filt) 


Le développement des termes en cosinus et en sinus permet d’obtenir : 



9 = — arctan 


2 SQ 



Régime 

transitoire! 


Régime permanent 


◄ ►.< 


x(t) 


temps 


Cas où A(t) est une fonction sinusoïdale 
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On remarque que dans tous les cas la solution homogène tend zéro lorsque t augmente, et 
la solution générale s’identifie alors avec la solution particulière : 

y(t) ^ yp(t) 

Pour cette raison la solution de l’équat-ion homogène est appelée solution transitoire tandis 
que la solution particulière est appelée solution permanente. 

A. 3. 4 Cas où A(t) est une fonction périodique du temps 

Introduction 

Nous avons étudié dans le paragraphe précédent la solution de l’équat-ion différentielle dans le 
cas harmonique, c’est-à-dire lorsque le second membre est une fonction sinusoïdale de la variable 
t. En considérant le cas de fonctions périodiques, nous procéderons à une généralisation du cas 
harmonique. 

Développement en série de Fourier d’une fonction périodique 

Les fonctions périodiques Une fonction f(t ) est dite périodique, de période T, si pour 
tout t, f(t + T) = f(t ) où T est une constante positive. La plus petite valeur non nulle de T 
est appelé la période de f(t). 

Exemples : 

1. La fonction sin(7ri) reprend les mêmes valeurs pour t = 2,4,6 puisque sin[7r(t + 2)] = 
sin[7r(f + 4)] = sin[7r(f + 6)] = sin(f). Cependant 2 est la période de sin(7rf). 

2. La période de sin(Uf) ou de cos (Qt) est 27 r/U. La période de sin (nQt) et de cos (nQt) est 
2n /nVL. 

3. Une constante a pour période n’importe quel nombre positif. 

4. Les figures ci-dessous représentent deux exemples de fonctions périodiques non sinusoï- 
dales. 



Exemples de fonctions périodiques 
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Définition des séries de Fourier Soit f(t) une fonction périodique de période T, c’est-dire 
telle que fit) = fit + T), la série de Fourier ou le développement de Fourier qui correspond 
f(t) est définie par : 

OO 

fit ) = d- a n cos(nOt) + b n sin (nflt) 

Z 

n= 1 

où Q = 2n/T, n est un entier positif, a n et b n sont les coefficients de Fourier. Ces coefficients 
sont définis par les expressions suivantes : 


& n 

bn 


— / f(t)cos(nflt)dt 

-L ./n 


- / f(t)sm(nQt)dt 


La valeur a 0 /2 représente la valeur moyenne de fit) sur la période T. Le terme sinusoïdal de 
pulsation f fi = Q . la plus faible est appelé fondamental tandis que les termes de pulsation 
Q n = nQ plus élevée sont appelés les harmoniques. 


Cas des fonctions paires et impaires 

— Une fonction est dite paire si f(—t ) = f(t). Dans le cas de développement en série de 
Fourier d’une fonction paire, il n’y aura que les termes en cosinus et parfois une constante 
qui représente la valeur moyenne. Il est aisé de montrer que : 


b n = 0 
4 

a n ~ Tp 

- Une fonction est dite impaire si f(—t ) 
sents dans un développement en série 
aisé de montrer que : 


( 2 r , \ ( 27mA 

^ f(t) cos ( - ~y~ J dt 

= —fft). Seuls les termes en sinus peuvent très pré- 
de Fourier d’une fonction impaire. Il est également 


ci n — 0 

7 4 A / 2nnt\ 

bn = T J Q sm \~Y~) dt 


Solution de l’équation différentielle La fonction A(t) étant périodique, de période T, 
son développement de Fourier s’écrit : 

OO 

A(t) = -^ + a n cos (nflt) + b n sin (nÜt) 

Zi 

n — 1 

L’équation différentielle s’écrit alors : 


ÿ + 25ÿ + cUoi/ 


a o 

~2 


OO 

+ a n cos (nQt) + b n sin (nQt) 

n = 1 
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La réponse permanente ( ou stationnaire ) qui s’identifie avec la solution particulière, pour 
t suffisamment élevé, peut alors être calculée pour chacune des composantes de l’excitation : 
oq/2, a.„cos(nOt), b n sin(nOf). On obtient alors par superposition : 


y(t) = 


a o 

2Ô| 


E 

n = 1 


a n cos (n n t + 0 n ) + b n sin (fl n t + 9 r 

- w §) 2 + « 2sj . 2 . 
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